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Final  Technical  Report 
(7/1/89  -  6/30/94) 

Abstract 

This  final  report  is  a  summary  of  all  of  the  research  performed  on  the  project  entitled 
Termolccvlar  Association  of  Ions  m  Gases  under  AFOSR  Grant  no.  AFOSR-89-0426  for  the 
period  7/1/89-6/30/94.  Theoretical  research  was  completed  and  published  on  the  following 
projects: 

(A)  Termolecular  Recombination 

A  +  B  +  U  —  AB  +  U 

(B)  Laser-Assisted  Electron-Excited  Atom  Collisions 

«”  +  X+  NKu  — » +  A *  +  MKu 

(C)  Electron-Excited  Atom  Collisions 

«”  +  A  -*  «“  +  A*,  A  =  B,  Be 

« 

(D)  Atom-Excited  Atom  Collisions 

A  +  B(n)  — *  A  +  B*  + 1~ 


A  total  of  23  publications  were  obtained  together  with  four  Ph.D  theses  during  the 
period  of  the  Grant.  Full  reference  to  this  work  is  provided.  In  addition,  new  theories  of 
Ion-Molecule  Collisions  (E)  and  Dissociative  Recombination  (F)  have  also  been  developed 
and  are  provided  in  this  report  as  well.  Full  details  of  the  theory  of  ion-molecule  collisions 
are  provided  in  an  accompanying  Organisation  Report  GIT-89-023. 
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X.  Introduction 

Ilia  Sad  report  summarises  all  of  tie  research  performed  oa  the  project  «a titled  TtrmoUcnUr  Asso¬ 
ciation  of  Ions  ta  Cum  dariag  the  period  of  Jeljr  1,  1989  through  Jaae  SO,  1994.  The  present  theoretical 
research  was  performed  aader  the  auspices  of  AFOSR  Great  ao.  AFOSS4MtU. 

The  objectives  of  the  research  program  was  to  formulate,  develop  aad  implemeat  new  theoretical  de¬ 
scriptions  of  various  atomic  aad  molecular  processes  of  importaace  ia  various  situations  of  interest  to  the 
Air  Force. 

XX.  Research  Completed 

Theoretical  research  eras  completed  oa  the  following  projects: 

(A)  Termdecular  Recombination 

A  +  B  +  M  ->AB  +  M 

(B)  Laser-Assisted  Electron- Excited  Atom  Coffiskms 

e*  +  A  +  Nku  -  e"  +  A*  +  M hu 

(C)  Atom-Excited  Atom  Collisions 

A  +  B(»)-*A  +  R+  +  e- 

(D)  Ion-Molecule  CoOirioas 

X*  +  AB  -*  {XAB+Y  -  products 

JT+  +  ABC  -  (XABC+y  -  products 

(E)  Electron-Excited  Atom  CoHiaioBa 


s'  +  A  -*  e“  +  A*,  A  =  H,St 

(F)  Empirical  and  Semiempirical  Representations  of  Ion- Atom  and  Atom- Atom  Interaction 
Potentials 

(G)  Electron-Ion  Dissociative  Recombination 


*-  +  AB+  ~>A  +  B 
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Each  of  the  Topics  ( A)-(D)  won  ik«  subjects  of  fou  Pb.D  thesis  performed  udci  the  auspices  of  AFOSR 
during  the  grost  period.  The  research  was  writtea  sp  osd  published  is  wiou  referred  jouub-  Reprists 
oad  Asses!  Reports  were  seat  routinely  to  AFOSR  eader  report  nos.  GIT-89-001  through  GIT-89-015.  Is 
odditios  to  this  prsssat  report  (so.  GIT-69'024),  the  reports  (GIT-88-016)  through  (GIT-88-023),  costsis 
o  copy  of  the  most  recast  Ph.D  thesis  oad  the  remaining  reprints  of  published  research  ore  attachments  to 
this  bask  final  report. 

Ill.  Published  Research 
A.  Publication  in  Reviewed  Journals: 

Preprints  of  the  following  research  were  submitted  to  the  AFOSR  under  the  GIT  Reports  specified 
below. 

1.  Maasky  E.  J.  and  Flannery  If.  R.,  The  Issue  of  Rests  Sot  Six*  »e'  +  F(li  -*  2e,2p) 

Collisions,  J.  Phys.  B:  At.  lfoL  Opt.  Phys.  28  601-7  (1980)  (Report  GIT-89-004). 

2.  Maisky  E.  I.  sad  Flannery  If.  R.,  Polarisation  Fractions  for  1ke2lP,SlP  *nd3lD  States  of  Helium, 
J.  Phys.  B:  At.  lfoL  Opt.  Phys.  38  3987-92  (1990)  (Report  GIT-89-005). 

3.  Maasky  E.  J.  and  Flannery  If.  R-,  The  JfuRschessel  Eikonol  Theory  of  Electron- Hydrogen  Collisions 
I.  Excitation  of  H(lt),  J.  Phys.  B:  At.  MoL  Opt.  Phys.  28  4549-72  (1990)  (Report  GIT-89-00fl). 

A  lfaasky  E.  J.  and  Flasasry  If.  R,  The  Multichannel  Bikonsl  Theory  of  Electron-Helium  Collisions 
I.  Excitation  of  J.  Phys.  B:  At.  lfoL  Opt.  Phys.  38  4573-4804  (1990)  (Report  GIT-89-007). 

5.  Smith  P.  H.  G.  and  Flannery  If.  R.,  Electron-Atom  Collision*  tu  a  Laoor  Field,  NncL  Instr.  Meth. 
Phys.  Res.  B  58/47  168-9  (1991)  (Report  GIT-89-011). 

6.  Smith  P.  H.  G.  and  Flannery  If.  R.,  Electron-Hydrogen  Collision t  ta  a  Loser  Field,  J.  Phys.  B:  At. 
MoL  Opt.  Phys.  24  L489-94  (1991)  (Report  GIT-89-010). 

7.  Maasky  E.  J.  and  Flannery  M.  R.,  Indirect  Coupling  Mechanisms  and  Stokes  Parameters  for  Electron- 
Atom  Scattering,  ].  Phys.  B:  At.  MoL  Opt.  Phys.  24  L551-8  (1991)  (Report  GrT-89-013). 

8.  Flannery  M.  R.,  Transport- Collision*!  Master  Egnations  for  TCrmolecnlar  Recombination  as  a  Func¬ 
tion  of  Gas  Density,  J.  Chem.  Phys.  95  8205-26  (1991)  (Report  GIT-89-012). 

9.  Smith,  P.  H.  G.  and  Flannery  M.  R.,  Electron-Hydrogen  Collisions  with  Dressed  Target 

and  Voiko*  Projectile  States  in  a  Laser  Field,  J.  Phys.  B:  At.  MoL  Opt.  Phys.  25  1021-49  (1992) 
(Report  GIT-89-009). 

10.  Maasky  E.  J.  sad  Flannery  M.  R.,  Eleetron-Meiaetable  Helium  Differential  and  Integral  Cross 
Sections,  2.  Phys.  B:  At.  MoL  Opt.  Phys.  25  1591-7  (1992)  (Report  GIT-89-016). 

11.  Flannery  M.  R.,  Termoleeular  Ion-Ion  Recombination,  Acta  Physics  Uaiversitatis  Comenianal  88 
11832  (1992)  (Report  GTT-89-017). 
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12.  Musky  E.  I-  ud  huwjr  M.  R-,  Empiric*!  and  Semiempirioal  /uterutim  Potentials  for  JUn 
Get- Ren  Gee  end  Mere  Gas-BaHie  Systems,  J.  Chtm.  Phys.  99  1902-77  (1993)  (Report  GIT-99-018). 

13.  BaSad  A.  ud  Flannery  M.  R~,  Angular-Momentum  Transfer  n  Colksiouel  Ionisation,  Phys.  Her. 
A  SO  429-39  (1994)  (Report  GIT-99-019). 


B.  Chapters  Published  is  Books  during  Grant  Period) 


14.  Flannery  M.  R.,  Recombination  Processes,  in  Molecular  Processes  is  Space  T.Wnianabe,  I.  Shima- 
mora,  M.  Skimin  aad  Y.  Itikava  (eds.),  Pleaam  Press  (1990)  (Report  GIT-99-002). 

15.  Maasky  E.  J.,  Electron  Collision  Cross  Sections  hoolwinf  Bodied  States,  ia  Nonequilibrium  Pro¬ 
cesses  m  PertUUp  Ionised  Geese,  NATO  ASI  series  B  320  349-55,  M.  CapiteQi  aad  J.  N.  Bardsley  (eds.), 
Pleaam  Press  (1990)  (Report  GIT-99-003). 

10.  Flannery  M.  R.,  Microecopic  end  Mecroecopic  Theorise  of  Termoleculer  Recombination  betvee a 
Atomic  lone,  in  Dissociates*  Recombination:  Theory,  Bepcriment  end  Applications,  NATO- ASI  series  B  913 
205-19,  B.  R.  Rowe  aad  J.  B.  A.  MitckeO  (eds.),  Pleaam  Press,  NY.  (1993)  (Report  GIT  89-020). 

17.  Flannery  M.  R.,  Electron-Ion  end  Ion-Ion  Recombination  Processes,  Advances  In  At.  MoL  Opt. 
Phys.  93  117-47,  Academic  Press  (1994)  (Report  GIT-89-021). 

18.  McDaniel  E.  W.  and  Maasky  E.  J.,  Guide  to  Bibliographies,  Books,  Unions  end  Compendia  of 
Beta  on  Atomic  Collisions,  Aduencee  ia  At  MoL  Opt  Pipe.  Special  Issue:  Cross  Section  Data  33  389-403, 
Academic  Press  (1994)  (Report  GIT-99-022). 

C.  Book  Barlows  Published  during  Grant  Pariodr 

e 

1.  Maasky  E.  J.,  Charge  Exchange  end  the  Theorg  of  Ion  Atom  Collisions  by  B.  H.  Bransden  ud  M. 
R.  C.  McDowell,  Oxford  U.  Press  (1992),  reviewed  in  Pkyrics  Today  40  124-25,  October  1993. 

D.  Personnel  Involved  during  the  Grant  Period) 

e  Graduate  Students  awarded  Ph  J)  degrees  during  Grant  Period 

1.  Dr.  M.  S.  Keenan,  Pk.D  thesis:  Termoleculer  Ion-Atom  Association,  awarded  3-17-90,  US  arisen. 

2.  Dr.  PkUUp  H.  G.  Smitk,  Pk.  thesis:  A  Semiclastical  Treatment  of  Laser  Assisted 
Collisions  in  s  Soft-photon  Week-field  Regime,  awarded  0-3-91,  non-US  arisen. 

3.  Dr.  A.  Hafisd,  Pk.D  thesis:  Angular  Momentum  Transfer  in  Electron-Atom  end  Atom-Atom  Colli¬ 
sions!  Ionisation,  awarded  8-15-91,  non-US  arisen. 

4.  Dr.  X.  Qi,  Ph.D  thesis:  Ion-Molecule  Spiraling  Collisions  end  Termoleculer  Recombination  (Report 
GIT-89-023),  awarded  9-16-94,  non-US  arisen. 


s  Senior  Research  Scientists 


1.  Dr.  E.  J.  Maasky  II,  Senior  Research  Scientist,  US  arisen. 
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(A)  Termolecular  Recombination 


IT.  Rmuch  Highlights  t 


A  +  B  +  M  —  AB  +  M 

Scientific  Objectives 

(1)  Develop  the  first  fundamental  coapteheaaive  microscopic  theory  of  this  simplest  three-body  chemical 
reaction. 

(2)  The  theory  to  serve  as  a  case  study  or  proto-type  tot  more  complex  three  body  transport-influenced 
reactions. 

(S)  The  rate  of  the  reaction  mast  be  famished  as  a  function  of  the  gas  (atomic  os  molecular)  species 
Jf  and  most  illustrate  the  non- linear  variation  of  the  rate  between  the  reaction  limited  regime  at  low  gas 
densities  and  the  transport  regime  in  the  Emit  of  high  gas  densities. 

Approach: 

(a)  Develop  sets  of  transport-coQisioaal  master  equations  which  govern  the  behavior  of  the  microscopic 
distribution  *(£,  E,  L)  at{A-B)  pairs  in  the  gas  if  over  then  internal  separation  R,  relative  energy  E  and 
relative  angular  momentum  L. 

(b)  The  set  of  equations  incorporate  a  blend  of  statistical  mechanics  and  theories  to  describe  the  trans¬ 
port  of  A  towards  B  through  the  gas  if ,  and  atomic  and  molecular  scattering  theory  to  describe  the  reactive 
and  non- reactive  collisions  between  the  pure  (A  -  B)  with  gas  spades  if. 

Accomplishments: 

A  major  paper  on  the  construction,  development  and  solution  of  the  required  equations  has  been  pub¬ 
lished  in  J.  Chem.  Phys.  05  (1891)  820&-26.  Publication  no.  5  of  Jill.  Reprints  sent  to  AFOSR  as  Report 
GIT-89-012.  Additional  papers  on  Recombination  which  have  been  published  ate  nos.  11, 16  and  IT  of  fill. 
The  latter  three  papers  have  eleo  been  seat  to  AFOSR  as  Reports  GIT-89-01T,  -820  and  -021,  respectively. 

Fundamental  Interest: 

This  development  has  served  as  a  textbook  study  of  the  most  basic  three  body  chemical  reaction  from  a 
microscopic  viewpoint.  It  provides  quantitive  and  physical  foundations  for  various  macroscopic  treatments. 
It  illustrates  how  reaction  and  traasort  are  coupled  in  transport-influenced  reactions. 

Applications: 

The  termolecular  proccse  is  key  to  the  basic  understanding  of : 

•  Rare  gas  *  halide  (exdplex)  lasers  and  rare  gas  exdmer  lasers, 

•  combustion  and  rocket  plumes, 

•  various  low  temperature  plasmas,  plasma  etching  and  plasma  decomposition  processes,  etc. 
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(B)  Laser* assisted  Electron-Hydrogen  Collisions 


e"  +  A  +  a hv  — • «“  +  A*  +  mhv 


Scientific  Objectives: 

(a)  To  develop  Ute  first  coapnhnsvi  theory  of  the  coBuioa  is  which  the  states  of  the  atom  A  an  drseeed 
(modified)  by  the  laser  aad  an  thea  closely  coupled  by  the  electro a-atom  iateractioa.  The  projectile  dectroa 
is  also  affected  by  the  field  of  the  laser. 

(b)  To  investigate  to  what  exteat  the  preseacs  of  the  laser  eahaacss  aad  caa  coatrol  the  degree  of  exdtatioa 
ia  electros  atom  coBmbobs. 

Approach: 

A  semklaaeical  Floqaet  approach  is  seed  to  fiad  the  dressed  states  of  the  atom  ia  the  laser  field. 
Volkov  states  an  the  states  of  the  dectroa  ia  the  laser  field.  A  semklaeskal  multichannel  dhonal  treatment, 
previoady  developed  for  electron-atom  coIBaoas,  thea  aeee  the  Floqaet  aad  Volkov  states  ia  a  doedy  coapled 
calcnlatioa  to  prodace  the  croee  sectioas  for  dectroa-atom  coHnioBS  ia  a  laser  fidd. 

Accomplishments: 

Sack  a  theory  was  developed  aad  applied  to  doctroa-hydrogea  colfaioM  ia  the  fidd  of  a  CO%  laser. 

e"  +  if(ls)  +  nhv  -*  e"  +  E (2e,  Ip)  +  mhv 

The  three  papers  (aos.  5,  6  aad  9  of  |XXX)  wen  published  aad  npriats  Mat  to  AFOSR  aa  Reports 
GIT-89-019,  -010  aad  -Oil,  respectively.  A  Ph.D.  thesis  entitled  A  SmiclassicaJ  Trt otmeni  of  Lcser-sssuted 
Collisions  in  e  Soft-pkoton  Wtmi-fisld  Regims  was  awarded  to  P.  H.  G.  Smith  ia  1991  under  the  supervision 
of  M.  R.  Flannery. 

(C)  Atom-Excited  Atom  Collisions 

Research  has  been  completed  oa  compatiag  tbs  cross  sectioas  lor  aagular  momeatum  changes, 

A  +  B(«f)-*  +  B+ +«'(«, O 

ia  heavy-particle  aad  dectroa-atom  («"  -  B)  coQisioas  when  the  target  atom  is  initiaDy  ia  highly  excited 
states  (n  >  10). 
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The  inuitk  ku  ilm  Utt  the  cross  sections  for  ni  -*  %t  cnllWwnl  transitions  iieniM  as  (  ia 

mcnind  until  ft  maximum  £ _ is  attained  After  which  the  cnee  sections  decreaee  precipitously.  TVie  effect 

eot  only  cae  be  — bp  *  qftftfttBJD  description  bat  a*  ahova  here  by  classical  scattering .  A  paper 
deecribiac  this  research  was  published  ia  Phys.  Rer.AU  429-39  (1294)  (ao.  12  of  Jill).  A  Ph.D.  then 
entitled  Anfuiar  Momentum  lVaa aftr  ia  Electron-Atom  sad  Atom-Atom  ColUnonal  Ionization  eras  awarded 
to  A.  Haffed  ia  1991  aader  the  supervision  of  M.  R.  Flannery.  The  above  paper  and  thesis  have  been  seat 
to  AFOSR  as  Report  GIT-89-019. 

(D)  Ion-Molecule  Collisions 

Ia  the  coatee  of  research  oa  recombiaatioa  it  was  realised  that  cross  sections  were  required  for  spiralhag 
collisions  betweea  ioas  aad  aeutral  molecules  with  permanent  dipole  moments  such  as  HCl  (linear)  aad 
NHf  (symmetric  top).  The  ion-molecule  interaction  is  ao  loafer  spherical  bat  is  orientation  depeadeat  eg. 
at  long  range  the  interaction  is, 

V(R,«)  =  ~  +  ^eos8  (3) 

where  a  is  the  polarisabiUty  of  the  molecule  with  permanent  dipole  moment  D,  and  #  is  the  angle  betweea 
the  R-axis  joining  the  ion  to  the  csnterwrf-maas  of  the  molecule. 

The  main  problem  in  these  ion-molecule  collisions  is  how  to  address  the  rotation  of  the  molecule  during 
the  course  of  the  coHirioa.  During  the  past  decade  there  have  beea  several  approaches  advocated. 

e  Semiclattical 

(a)  The  Loch-in  Dipole  Approximation 

(b)  The  Frosea  Rotor  Approximation 

(c)  The  Average  Dipole  Orientation  Theory  (ADO) 

(d)  The  Free  Energy  Average  Angular  Motion  Theory 

(e)  The  Adiabatic  Invariance  Theory  (ADIA) 

s  Quantum  Mechanical 

(a)  The  Pert  or  bed  Rotational  State  Theory 

(b)  The  Adiabatic  Capture  and  Centifrugal  Sudden  Approximation  (ACCSA) 

One  of  the  results  of  the  present  research  has  been  the  modification  of  the  adiabatic  invariance  theory, 
to  account  for  the  coupling  betweea  the  internal  angular  momentum  of  the  target  molecule  aad  the  orbital 
angular  momentum  of  the  projectile  about  the  target,  in  Ac  computation  of  the  rate  k  of  spiralling  collisions 
ia  systems  as  -  BCl,  a  rate  analogous  to  the  Langevin  temperature- independent  rate, 

g 
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*X  =  2.M.10-*(£)l/*em*«-1  (4) 

for  spiraffiag  cflBhwn  udtr  Ik*  polarisation  attraction  ilou  ie.  Ik*  int  Una  of  (3).  Is  (4)  a  it  Ik* 
polarisabOity  is  mrnht  of  A*  and  tk«  ion-molecule  ndaced  aaat  if  it  is  am*.  Raccat  experiment  (D.  C. 
Clary,  D.  Smilk  mad  N.  O.  Adana  Cktm.  Phys.  Lett,  lit  (IMS)  320)  kaa  akowa  tkal  tkcrc  it  a  large 
difference  between  k  (T)  which  nmt  witk  iaotkemal  Icnptralan  T  amd  Ik*  Laagevia  rat*  hi  vkick  it 
temperature  independent.  Cakalalioa  of  Ike  rat*  k  it  quit*  complicated  even  Cm  systems  at  Jfe+  -  SCI 
aad  Bf  -  SCI  by  nk- initio  qaaalan  nackaaical  lh*ori*t.  A  Pk.D  thesis,  eatitlad  Jan-Moltnlt  SjnrnBnf 
Collision*  and  TtrmoUenUr  Rtcomkinoiion,  vat  awarded  to  X.  Qi  ia  1004,  aad  kaa  kata  aaat  to  AFOSR  at 
Report  GIT-S0-023. 


(B)  Electron-Excited  Atom  Collisions 

Ia  tk*  developmeat  of  semirlassifal  tkaoriaa  for  coffitiot  procatttt,  «"  +  A?  -» «"  +  AJ,  wkere  i  aad  j 
ate  botk  excited  or  metaatable  electronic  atatea  of  tke  atom  A ,  tke  practical  implcmeatatioa  of  fully  qaaatal 
theories  evea  if  desired  is  unfeasible  with  modem  supercomputers.  Somekow  tke  pkytict  emeatial  ia  mllitioat 
witk  excited  atom*  it  iadBdeatfjr  daacribed  by  modera  qaaatal  metkoda  ia  tk*  teat*  tkat  ia  order  to  obtaia 
results  considered  to  be  accarate,  toar-de-ibre*  eraelatioa  of  a  wkolt  koet  of  email  terma  ia  reqeired.  Ia  feOy 
qaaatal  cakalatioaa,  large  amoaata  of  compatcr  time  are  ipeat  evaluating  termi  wkick  ehimately  provide 
insignificant  coatribntioa  to  tke  cram  section.  Ia  felly  qaaatal  approackaa  tkere  ia  ao  a- priori  mrtkod  of 
tk*  regioa  wkick  effectively  controls  tke  cram  aectioa.  Semklaaaical  teckaiqam  tkerefore  kaa*  beta 
developed  aad  aaed  to  track  tke  emeatial  pkyaica  mack  more  efficeatly  aad  effectively.  Cam  atadim  for  tke 
target*  hydrogen  aad  kdiam  kave  beea  carried  oat. 

Tkia  work  kaa  pabfiaked  (ao*.  1-4,  7,  10  aad  15  ia  Jill)  aad  reprint*  km  been  aeat  to  AFOSR  m 
Report*  GIT-8»>004,  -005,  -006,  -007,  -013,  -016  aad  -003,  respectively. 

(F)  Empirical  and  Semiempirical  Representation*  of  Ion- Atom 
and  Atom-Atom  Interaction  Potentials 

Tke  Ttag-Toeaaiee  (TT)  mmi-empirical  modd  potential*  for  ion-atom  eyatema  i*  applied  to  tke  rare 
gaa-halide  negative  ion  exripkrm.  Tke  coefBdent*  defining  tke  repulsive  Born- Meyer  term  ia  tke  TT  semi- 
empirical  potentials  are  determined  from  the  equilibrium  bond  length,  and  dissociation  energy,  D,,  taken 
from  nh-initio  calculations  aad  from  transport  studies  at  them  molecular  ions.  The  damped  dispersion  aad 
induction  energy  terms  in  the  TT  potentials  are  obtained  from  coupled  Hartreo-Fock  calculations  for  the 
neutral  rare  gas  atoms  aad  F~,  Cl~  ions.  The  multipole  polarisabifitim  for  the  heavier  halogen  atomic  nega¬ 
tive  ioas  are  estimated  from  a  knowledge  of  polarisability  ratio*  across  isoelectronie  sequences.  The  resultant 
semi-empirical  ionic  potentials  axe  compared  to  available  at- initio  calculations  aad  the  results  of  inversion 
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of  tnuport  tWy.  To  facfliUte  the  conpuna  of  tke  (spam)  ak-mttio  data  vttk  Ike  semi-empirical  po- 
Uatiah,  a  uph  |t(ia|  pncadan  k  prsa—tsd  (m  determining  empirical  pot— tiak  ioi  diatomic  molecules 
from  a  ml  of  Am  coutniat  equations.  Tko  fttting  procadare  k  applied  to  a  total  of  33  nit  gas  — rimoto 
aad  nit  gaa-kabdt  — dpi— m  (botk  n— teal  and  ink)  of  iaUroot  to  a  variety  of  applications  in  gassoaa 
discharges  aad  exdaet  1 — rn  A  Vterm  aopao— tarin  of  tko  empirical  pot— tiak  g— orated  k  accomplished 
with  tko  an  of  a  data  eet  which  iadade  tke  ‘geometric’  panmeteie  (Ao,  A«,  D,}  aad  the  additional 

fMurameteae  (a*,  I.P.,  EJL }  needed  for  tke  disperrion  aad  iadnctioa  energy  terms.  A  no— 1  foatare  of  tke 
empirical  procedare  k  tke  form  elation  of  tke  coaeteaiat  equations  at  two  and— r  dkplacem— to  (1  constraint 
at  Rq  vkcrtia  tke  potential  paeem  tkroagk  aero,  and  3  constrain  to  at  A*,  tke  equilibrium  separation)  vkkk 
yields  an  accarate  ft  to  available  ak-tasrio  data  aad  greatly  extends  tke  mage  of  in  ten-dear  separations 
R  for  vkkk  an  accarate  piecevke  analytical  empirical  potential  can  be  generated.  Tb  teat  tke  relative 
importance  of  tke  difereat  terms  in  tke  ftted  3* term  empirical  teptee— tationa,  tke  daaairal  orbiting  croea 
section  Q«ru<(£)  is  competed  aaing  tke  faO  empirical  potential  aad  compand  against  tke  standard  Laagevin 
orbiting  cron  section  Qp*.(E)  tot  a  pare  polarisation  interaction. 

Tkk  work  has  bo—  pablkkod  in  J.  Chem.  Pkya.  (no.  11  of  |ZZX)  and  reprints  have  be—  ant  to 
AFOSR  as  Report  GIT-88-018. 


(G)  Dissociative  Recombination 

Ahhoagk  aO  tke  modem  qaantnm  (scattering  aad  chemistry)  technology  has  be—  broagbt  to  bear  — 
dissociative  recombination, 


«“  +  AB*  -*  A  +  B  (2) 

for  tke  simplex  diatomica,  then  remains  several  incoukteadee  with  observational  data 

A  asv  dam  of  diaaodative  recombination  (DR)  k  emerging.  In  contrast  to  normal  DR,  cknracterised 
by  rates  apx(T)  —  3  •  10~,(300/T)1^*  m  for  Of  (if  =  0)  and  NO*  vitk  large  dkaodetive  —  ergiee  Df  ~  7 
eV  — d  11  eV,  a  dam  characterised  by  —per  rates  of  3  •  10~a  are  being  discovered. 

Therefore  a  aaad  to  investigate  dimodativs  recombination  for  varioas  complex  systems  by  tke  develop¬ 
ment  of  physical  theories  vkkk  viD  furnish  insight  into  the  various  mechanisms  k  varr anted  by  tke  ongoing 
inconsiet— dm  vitk  experimental  data.  Tke  sk-taftto  calculations  are  having  n  tough  enough  time  even  for 
the  simpler  specks  Bf,  Nf,  Of,  etc.,  and  are  totally  impractical  tot  more  complicated  systems  of  interest. 

A  review  of  electron- ion  — d  ion-ion  recombination  processes  has  recently  be—  published  (no.  17  of 
|ZZZ).  A  reprint  of  this  review  has  be—  a— t  to  AFOSR  as  Report  GIT-88-021.  Among  the  topics  covered 
are  the  new  clam  of  dissociative  recombination  reactions  mentioned  above.  Progress  made  towards  the 
developm— t  of  a  new  theory  k  piss— ted  in  of  tkk  report. 
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f .  Ckidtil  Path  ThMrjr  of  Direct  Electron- Ion  Dissociative  Recombination 

and  Associative  Ionisation 


T.l  Background 

In  1950  Bates  [1]  postulated  that,  dissociative  recombination  (DR)  Cor  diatomic  ions  can  occur  via  a 
crossing  at  Rx  between  the  bound  and  repulsive  potential  energy  curies  V+(R)  and  for  AB*  and 

AB*‘,  respectively.  Here,  DR  involves  the  two-stage  sequence, 

ht  V, 

«-  +  AR+(ei)  ee  (AB**),  —  A  +  B*  —  A+B+fc* 

The  first  stage  is  dielectromc  capture  whereby  the  free  electron  at  energy  «  =  V**(R)  -  V*(R)  excites 
an  electron  of  the  diatomic  ion  AB*  with  internal  separation  R  and  is  then  resonantly  captured  by  the 
ion  at  rate  h,  to  form  a  repulsive  state  d  of  the  doubly  excited  molecule  AB**,  which  in  turn  can  either 
autoionise  at  probability  frequency  v«,  or  else  in  the  second  stage  predisaodate  into  various  channels  at 
probability  frequency  v,.  This  competition  continues  until  the  (electronically  excited)  neutral  fragments 
accelerate  past  the  crossing  at  Rx-  Beyond  Rx  the  increasing  energy  of  relative  separation  has  reduced  the 
total  elsrtronir  energy  to  such  an  *  that  »*«**«t<*fliiit is  essentially  prsrlmisd  *****  neutralisation 
is  then  rendered  permanent  past  the  stabilisation  point  Rx-  Bates'  interpretation  has  remained  intact 
and  robust  in  the  current  light  of  s5-ta«h's  quantum  chemistry  and  quaatal  scattering  calculations  for  the 
simple  diatomics  (Of ,  Nf ,  Naf ,  etc.).  Observation  of  emitted  radiation  h»  yields  information  on  the  excited 
products.  Mechanism  (1)  is  termed  the  direct  process. 

In  1998  Bardsley  [2]  pointed  out  the  possiblity  that  a  three  stage  sequence, 

e*  +  AB+(ui)  -*  [AB+(v/) -  e-J„  -*  (AB**),  -» A  +  B*  (2) 

the  so-called  sndtrset  process  might  contribute.  Here  the  accelerating  electron  loses  energy  by  vibrational 
excitation  (ej  -*  ey )  of  the  ion  and  is  then  resonantly  captured  into  a  Rydberg  orbital  of  the  bound  molecule 
AB*  which  then  interacts  one  way  (via  configuration  mixing)  with  the  doubly  excited  repulsive  molecule 
AB**.  The  capture  initially  proceeds  via  a  small  effect  -  vibronic  coupling  (the  matrix  element  of  the  nuclear 
kinetic  energy)-  induced  by  the  breakdown  of  the  Born-Oppenheimer  approximation  •  at  certain  resonance 
energies  <«  =  B(e/)  —  B(v?)  and,  in  the  absence  of  the  direct  channel  (1),  would  therefore  be  manifest  by 
a  series  of  characteristic  very  narrow  Lo rents  profiles  in  the  cross  ssetion.  Uncoupled  from  (1)  the  indirect 
process  would  augment  the  rate.  Vibronic  capture  proceeds  more  easily  when  */  =  e<  +  1  so  that  Rydberg 
states  with  n  «  7-  9  would  be  involved  (for  Hf(a  =  0))  so  that  the  resulting  longer  periods  of  the  Rydberg 
electron  would  permit  changes  in  nuclear  motion  to  compete  with  the  electronic  dissociation.  Recombination 
then  proceeds  as  in  the  second  stage  of  (1)  ie.  by  electronic  coupl  sg  to  the  dissociative  state  d  a*  the  crossing 


point.  Giusti  [3]  has  ptmJiid  *  unified  account  of  th«  direct  and  indirect  proceeeee. 

O'Malley  [4]  later  noted  that  the  proceee, 

*c  Vi 

•~+AB+{*)  ~  (AB")<  — *  A  +  B*  (3) 

v, 

V*4  T1  VM 

[AB+(v)-4-  L 

proceed*  via  the  firet  (dielectronic  capture)  stage  of  (1)  followed  by  a  two-way  electronic  transition*  with 
frequency  and  v*a  between  the  d  and  n  state*.  All  (n,  *)  Rydberg  states  can  be  populated,  particularly 
thoee  in  low  n  and  high  •  once  the  electronic  d— n  interaction  varies  as  »'u  with  broad  structure-  Although 
the  dissociation  process  proceeds  here  via  a  second  older  effect  (i/*»  and  *»*)  the  electronic  coupling  may 
dominate  the  indirect  vibronic  capture  and  will  interupt  the  recombination  in  contrast  to  (2a)  which  as 
written  in  the  one-way  direction  foods  the  recombination.  Such  dip-structure  has  been  observed.  Guberman 
and  Giusti-Susor  [5]  have  assessed  the  effect  of  each  contribution  of  (1),  (2a)  and  (2b)  to  the  resonance  shape 
and  integral  cross  section. 

There  exist  two  a 3-mttto  quanta!  treatments  for  dissociative  recombination  -  one  (CM)  based  on  con¬ 
figuration  mixing  [2-4,6]  ,  and  the  other  (MQDT)  based  on  multichannel  quantum  defect  theory  [3].  CM 
has  been  applied  [6-9]  to  ;  MQDT  has  been  applied  [16-12]  to  H},  ,  to  NO*  [13],  to  Of  [6] ,  to  Nf 
[14]  and  to  CM*  [15].  These  a>-«attie  treatments  [10-12,15]  have  shown  that  the  indirect  process  interferes 
destructively  with  the  direct  process.  The  cross  sections  exhibit  an  overall  preponderance  of  destructive 
interference  via  a  series  of  dips  foiling  below  the  B~x^  continuous  variation  for  the  direct  process.  Here 
the  direct  — ♦  indirect  coupling  interrupts  the  recombination  at  specified  energies.  The  quantal  theory  for 
diatomic  ions  is  considered  to  be  essentially  complete  for  cases  involving  favorable  crossings  between  the  ion 
and  doubly  excited  neutral  states.  The  EtH*  case  which  does  not  involve  curve  crossing  has  been  treated 
recently  by  Guberman  [16]  and  by  Sarpal  si  s4  [17].  Although  the  crossing  at  ~  8  eV  for  Mf  has  been 
explored  by  Kulaader  and  Orel  [16],  the  large  cross  section  behavior  [19]  at  low  energies  (where  there  is  ae 
curve  crossing)  remains  unresolved. 

Chris  Bottchsr  [6]  in  1976  developed  a  rather  nice  semiclassicai  theory  of  dissociative  recombination. 
The  theory  is  smiirlassical  in  that  JWKB  wavefunctions  for  nuclear  motion  were  used  in  this  stationary 
state  quantum,  treatment.  Miller  [20]  has  examined  associative  ionnisation  •  the  inverse  of  (1)  •  within  a 
classical,  semiclaasical  and  quantal  framework.  In  this  paper  a  different  scmiclassical  theory  of  dissociative 
recombination  is  developed  from  a  time  dependent  theory  based  on  a  classical  trajectory  R  =  R(t)  for  the 
relative  motion  of  the  dissociating  neutral  heavy  particles. 
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A.  Motivation!  Why  Another  Theory  ? 

When  the  two  diabetic  potential  energy  curve*  V*(R)  and  t^(Jf)  cnee  at  Rx,  the  rate  for  dieeodative 
recombination  given  by  a  lint-order  treatment  [2,21]  is, 

— «  ■  jz£ivn(G)  [T^wiii^x]  { m 

where  l ^»(H)  is  the  electronic  energy,  (dj(JZ,r)  |  Hti(Rx)  |  +,(R,r))r  coupling  the  electronic  wavefunction 
de  for  the  doubly  excited  repulsive  (AR)?*  intermediate  state  to  the  electronic  wavefunction  for  the 
scattering  system  e  —  AB+.  The  molecular  functions  ds,<  arc  diabetic  in  that  they  an  net  pun  eigenstates 
of  the  full  electronic  fixed-nuclei  Hamiltoninn  Sti(R^r),  and  the  continuum  functions  an  energy-normalised 
(d«  I  d«*)r  *  ~  fl)  with  consequent  unit  density  f(t)  of  states  for  the  scattered  electron  of  energy  e.  The 

bound  vibrational  wavefunction  for  the  original  ion  state  is  j+(R). 

Although  all  the  modem  quantum  (scattering  and  chemistry)  technology  has  been  brought  to  bear 
for  the  simpler  diatomics  in  a  way  considered  in  general  correct,  simple  exprsas ions  as  (4)  an  invaluable 
in  that  they  reveal  tremendous  insight  into  the  mssntisi  physics,  general  characteristics  end  workings  of 
recombination  for  various  distinct  systems.  Bates  [21]  has  demonstrated  how  sxpreasious  as  (4)  may  be 
utilised  very  effectively  to  promote  new  insight  for  recombination  involving  more  complex  ions. 

Since  Bates’  deductions,  particularly  his  discovery  of  the  new  dees  of  super  dissociative  recombination, 
characterised  by  rates  as  large  as  2  •  10-*cm*s~l,  in  contrast  to  2  •  10~Tcm**_1  for  normal  DR,  rely  on 
.the  gradient  JVtfdR  of  the  doubly  excited  repulsive  state,  it  is  now  of  interest  to  sec  whether  (a)  simple 
expressions  as  (4)  but  more  accurate  can  be  derived  without  recourse  to  the  foil  numerical  approach  of  at- 
initio  treatments,  and  (b)  if  the  shape  of  the  ion  potential  V+  enters  as  directly  as  does  Vd  in  (4)  over  that 
implicitly  contained  in  ^(Rx)>  Apart  from  the  first-order  treatment  (4)  and  the  foil  quantal  treatments 
(CM,  MQDT)  above,  there  appears  to  be  no  intermediate  simplified  description.  There  are  also  no  analytical 
results  such  as  (4),  appropriate  to  the  case  whan  V+  and  Vt  do  not  cross. 

The  purpose  of  this  paper  is  to  provide  such  a  method  wherein  a  two  state  semiclassical  treatment  of 
the  direct  process  (1)  will  effectively  yield  an  expression  for  the  rate  in  a  form  similar  to  but  mote  accurate 
than  (4). 

V.2  Classical  Path" Theory 
A.  Basic  Equations 


In  the  field  of  the  ion  AB+,  the  recombining  electron  is  captured  at  time  t  =  t*  into  an  electronically 
doubly-excited  molecular  state  ds(r,  R)  whose  variation  with  R  is  provided  by  the  classical  trajectory  A(t ). 
Competition  between  autoionisation  to  the  combined  state  di(r,  JR)  of  the  emitted  electron  of  energy  *  and  the 


diatomic  ion  AB+{n)  continues  a*  the  nuclei  move  inward  or  outward  along  the  trajectory  R(t) 
with  the  potential  V^Jt).  The  system  wavefunction  satisfies. 


=  iA^M  (5) 

where  the  time  dependence  in  the  electronic  Hamiltonian  H^r,  R(t))  —  H~(h/2ii^f)Vjf  for  firoeen  nuclei  in 
terms  of  the  total  system  Hamiltonian  H,  is  generated  by  the  clsssical  trajectory  R(t)  tot  atomic  dissociation. 
Expand 


*««>■  £  e,(t)d.(r, R)omp [-i jT w^t))  eft]  («) 

in  terms  of  some  basis  set  dy  where  Wjj  are  diagonal  elements  of, 

i M*W)  =  {*(?,  *)!*.»  I  Mr,  J*)),  (7) 

Insert  (8)  into  (5),  project  onto  the  state  <  di  |  and  use  (^i  |  Ht)f  =  &  •  (d*  I  Vj,dj)  so  as  to  obtain  the 
set  of  coupled  equations, 

,  >  MO  ]  •»  (-ii,'  Mw>  -  Wtfwj  -)  (•) 

i+i 

tat  the  transition  amplitudes  e<  in  this  classical  path  representation.  When  dy  are  identified  with  the  adiabatic 
(molecular)  functions  xj  which  satisfy, 


Stouif, £)  =  [«*..(*)  +  V{f,  *)]  xj(r,  *)  «  *y(*)xy(f,  8)  (9) 

then  Wy(H)  =  Ej(R)6ij  in  (8),  curves  W«(.R)  and  Wjj(R)  of  the  same  symmetry  do  not  crocs  and  the 
i  —  j  electronic  transition  occurs  via  the  uon-adiebatic  dynamic  coupling  elements  (x»  I  V^xy)*  When 
dy  are  taken  as  (atomic)  eigenfunctions  of  the  internal  electronic  Hamiltonian  Hi*t  at  infinite  separation 
R,  then  these  dynamic  coupling  terms  vanish  and  the  transition  occurs  via  the  potential  coupling  terms 
Wy  s  (d<  |  V{ft  R)  |  dy  )•  When  dy  are  taken  eigenfunctions  of  some  Hamiltonian  Ho  intermediate  between 
H,i  and  Hi*t,  but  ao  chosen  as  to  satisfy, 

(d»  1 7<«t  |  di)y  *  A  »> 

*-<*|Vjildy>r~0  1  ’ 

where  7f«j  =  7<o  +  H!  with  H'  <  7<o,  then  dy  form  a  diabetic  basis  set  since  the  dynamic  couplings  are 
small.  A  complete  atomic  basis  set  di^a.rs.f*  -»  oo)  =  ds(r)ds(r)  centered  on  each  nuclei  A  and  B  (ie. 
Ho  =  Hi* t,  H'  =  V)  constitutes  the  amplest  diabetic  baas,  eigenfunctions  of  Hi*t  of  (9).  Application  of  (8) 
to  the  electronic  continuum  is  obtained  by  — p—«<«"g, 
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9(f,t)  =  c*(t)de(r.X)««P  [-J  jf  *„(t)  *]  +  c.(t)d.(r,X)e*p  *+(M)  dt]  *  (11a) 

in  terms  of  the  diabetic  orthogonal  bans  set  da(ri  R)  of  energy  normalised  functions  (d«  |  d«*)*  =  <(«  —  <') 
for  the  continuum  electron-ion  system  e~  -  vtB+ffl),  and  of  ds  for  the  molecular  system  AB**(R)  in  a 
vibrational  continuum.  The  energies  Wjj  an  the  electronic  potential  energy  surfaces, 


*+(«,*)  =  (d.  |  W.,  |  d^)  =  V+(H(t))  + 1  (Ilk) 

and, 


*e(i)  =  <d-!W.i|dc>  =  Vi(i)  (11c) 

within  which  nuclear  motion  proceeds.  Assume  now  that  all  continuum  states  are  all  uncoupled  ie.  (d«  | 
I  d«*)  =  d(<  -  «*).  In  this  diabetic  representation  (live),  then  (8)  yields, 


and, 


t  ft  ds(t)  =  j  Ku(<)c,(t)  *xprr(r,t)  dt 


(12a) 


*  ft  i,(t)  *  VJ,(t)c*(t)  asp  [— vy(e;  t)J . 


The  bound-continuum  electronic  coupling  matrix  elements  an, 


(12*) 


V*(t)  =  (d<  I  *«(r,R(t)}  |  d.(r,X))r,.  *  K,(t)  (12c) 

when  the  integration  is  over  the  electronic  coordinates  food  the  direction  f  of  the  ejected  electron.  The 
phase  is, 


+  *'  (12 d) 

The  above  clasrical  path  formulation  (12)  is  now  applied  to  the  forward  direction  of, 

At 

A  +  X*  AX++e"  (13) 

DR 

since  the  relative  A  -  X*  motion  can  be  described  by  the  classical  path  X(t)  and  the  electronic  motion  by 
quantum  mechanics.  The  strategy  is  therefore  to  solve  (12a)  for  es(t)  subject  to  c*(tx)  =  1  ie.  the  reaction 
does  not  begin  to  occur  until  the  A  -  X*  separation  is  Rx  at  t*.  Solution  of  (12b)  then  provides  the 
probability  for  electron  ejection  within  the  reaction  sone  R  <  Rx  ie.  for  the  decay  probability  of  a  discrete 


IS 


electronic  state  AB**  moving  via  the  classical  path  R(t)  through  as  embedded  electronic  continuum.  The 
crew  section  etju  for  the  forward  channel  (associative  ionisation)  of  (IS)  can  then  be  obtained.  Detailed 
balance  then  provides  the  cross  section  *dh  tot  the  reverse  channel  (dissociative  recombination)  of  (13)  ie. 
for  the  transition  involving  molecular  states  initially  in  an  electronic  continuum  and  finally  in  a  vibrational 
continuum. 

B.  Formal  Solutions 

With  (12b)  inserted,  the  formal  solution  of  (12a),  subject  to  the  initial  condition  e(t,tx)  =  0, 

-2#fiacs(t)  *  jf  dt'jf  d«r(«,  t;  «*F  *  W«i  *)  ~  7(<i  O)  (14) 

at  time  t.  This  depends  on  the  previous  history  of  the  system  between  tx  and  t  via  the  non-local  interaction, 

r(c,t;0  =  2»K«.(t)t^(0  (IS) 

and  on  the  phase  difference, 

y(e;t)-y(«;0  =  jjft[K.-(V++«)]  *  (16) 

at  different  times.  This  difference  f»"  be  expanded  to  yield, 


*«*>• -tfco  -  (£).« -«')+ i  (£)/* -*r+..  • 

=  j{IW)-v*(0]  -«}(i-t')+r<»-«'>  un 

=jW-<i(i-n+ni-fl 

where  F  is  (-independent,  being  a  function  only  of  the  difference  (t  —  If),  and  where, 

W(<)  =  Vt(t)  -  V+(t)  (18) 

is  the  energy  for  vertical  transitions  at  time  t.  Equation  (14)  with  (17)  therefore  reduces  to, 

-2*h*cs(t)  =  jf  eepiF(t-t')  dt' jf°  de  {r(c,  t;  1f)e4(tf)  exp  £  [W(t)  -  e]  (t  - 1')}  (18) 

which  is  in  a  form  suitable  for  further  approximation.  The  following  analysis  is  valid  for  the  case  V4(t)  > 
V+(i)  for  R  <  Rx  appropriate  to  curv*croesiag  between  the  ion  and  neutral  states  at  Rx  ie.  for  W(t)  >  0. 
The  separation  Rx  can  however  tend  to  infinity  so  that  VtfR)  >  V+(R)  everywhere. 
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C.  Local  Approximation 

The  «-  integration  ia  (19)  involves  the  product  of  T  •  e*  which  is  bow  assumed  to  vary  (lowly  with  «, 
(i*.  r  has  a  large  width  A  ia  it*  variation  with  c),  and  aa  exponential  which  oerillataa  over  c  with  period 
6  s  2wK/(t—f).  The  s-iategral  is  then  negligible  within  the  (t-t')-raags  which  satisfies  f  <  A  is.  for  t  >  f. 
The  remaining  range  0  <  (t  -  tf)  <  2wA/ A  will  therefore  provide  the  main  contribution  to  the  (-integral. 
Hence  e4(t)  of  (19)  has  only  a  short  memory  of  the  previous  values  of  c*  between  U  and  t  and  depends  only 
on  t'  immediately  before  t.  Thus, 

— 2wft*e 4(t)  m  e4(t)  £*  T((,  t)  dr  £  sap  [  j0*(0  -  «K]  *r  (20) 

for  all  t.  As  t  — » oo,  or  more  precisely  for  t  >  2*A/A  when  the  width  A/t  <  A,  then 

/'  -»  [k<-  — H  ■ * ■  ■ *  H-  -  ■ ■> +  ■  <JI> 

The  «—  integration  therefore  yields, 

*s(i)«-^t)[ir(«, «)+**,(«)]  (22) 

where  the  energy  width, 

r(*,t)  *  2v  |V*(t)|» ,  e  =  HT(i)  (23) 

is  due  to  coupling  between  de  and  one  continuum  state  d<  with  energy  H(t)  and  where, 


is  the  second-order  energy  shift  in  state  dr  due  to  the  coupling  with  all  other  continuum  states  d«(r.  R)  with 
energy  c  ^  W(t);  since  the  principal  value  part  V  involves  the  contribution  from  states  just  below  W(i)  to 
be  balanced  by  states  just  above  W{t).  The  solution  of  (22)  is, 

ce(t)  =  esp  [~L  £  T(t?)  «*]  esp  [-j£  SEtf)  *']  (25) 

The  probability  |c4|*  for  remaining  on  the  V4  curve  at  time  t  is, 

P4(<)  =  «ap[-i£r(0*']  Se.p[-^£|Vi.(H(<))|,H(<)df]  *  e.p  [-£  *.(<)*]  (26) 

which  is  the  probability  for  survival  at  time  t  against  ««^uTmiUttd  autoionisation  at  frequency  between 
tx  and  t.  The  effect  of  the  continuum  state  d«(r,  R)  on  the  dissociating  system  is  therefore  realised  by 
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assigning  the  compld  nogr  A*  (A)  *  Vt{  A)  4-  IBt(R)  -  }r(A)  to  state  4d  whose  occupation  probability 
then  decays  with  time  as  Pt(t).  Tht  solution  of  (12b)  ia  therefore, 

»M«,0=  f  C4(0Vi(0*vI-n(r, t'))A'  (27o) 

with  the  phone, 

(2Tb) 

expressed  in  terms  of  the  positive  energy 

W{t)  =  V^t)  -  V+(t)  (27e) 

When  the  local  eolation  (25)  for  c*(<)  is  inserted  into  (27a),  the  amplitude  c,(t)  can  therefore  be 
determined  by  direct  numerical  integration.  Initially  at  time  tx,  the  nuclei  A—B*  enter  the  (autoionieation) 
reaction  aone  at  Rx-  They  continue  their  relative  motion  inward  until  they  reverse  their  motion  at  A,,  the 
distance  of  closest  approach,  at  time  i«,  move  outward  and  then  leave  the  reaction  none  with  a  diminished 
amplitude  at  time  tj+Ty  with  eeparation  Rx-  On  taking  tx  =  0  and  t  =  rx  in  (27a),  the  amplitude  c,(t) 
for  electron  emission  after  the  time  interval  ry  following  some  rearrangement  is  given  by, 

t  h  c(«)  =  jf  es(t)VX(t)e*p[-n(ejt)]  «*  + [sap (t)VX(t)  sap  [+ry(<;t)]  dt  (28a) 
where, 

cs(t)  as  snp  J  T(<)  dtj  (286) 

is  the  amplitude  for  survival  of  travel  from  Rx  to  A(t)  on  the  inward  leg  (t  <  t(),  and, 

*s(t)  =  Ci[rx  -t)  =  e*p  [“Jf  W  dtj  esp  jf  T(t)  dtj  (28c) 

ia  the  amplitude  for  survival  during  the  sequence  Ax  ^A,"*  A(rx  —  t)  »  R(t)  on  the  outward  leg  and, 

’k’*)s=£jTiw,(<,",)  *s2yM  (2M) 

is  the  total  phase  accumulated  within  the  full  reaction  tone.  Note  that  e s(t)s*(t)  =  cs(rx)  =  cj(t.)  is  the 
survival  amplitude  for  the  round  trip  Ax  -*  R*  — 1 >  Rx-  The  above  form  (28a)  is  useful  in  that  it  tracks  the 
increasing  accumulation  of  phase  y(t)  which  in  (27b)  is  measured  from  the  time  <j  of  entry  to  the  reaction 
sane.  An  equivalent  but  symmetrical  form  of  (28a)  follows  by  taking  i  =  0  at  A,  such  that  <x  =  — t,  and 
t  =  Tx  =  +t«  in  (27a).  The  result  is, 


»fie(*)  =  e*p[-»A(t,)]jf  Ku(*)Mf)«*p{+«A(t)}  +  *(<) e»p{-»A(<)}l  A  (29a) 

when  the  new  phase, 


[*(*)- e]  A 


(29*) 


ia  sera  at  ft,  s  JZ(t  =  0).  Since  TT(t)  is  even,  the  old  and  new  phases  are  than  related  by, 


A(i.)-A(t), 
A(i.)  +  A(t), 


*<*. 

*>*. 


(29c) 


The  first  and  second  terms  in  each  of  (28a)  and  (29a)  are  the  respective  contributions  to  the  amplitude 
c(«)  from  the  inward  and  outward  legs  of  the  classical  trajectory  R (<).  Bach  form  is  useful  for  correspondence 


with  the  stationary-state  quanta!  transition  amplitudes  when  JWKB  bound  vibrational  wavsfunctions  are 
taken  relative  to  the  left-hand  or  right-hand  turning  points,  respectively  (cfc  |VZ). 

The  banc  classical  path  theory  is  represented  by  either  the  amplitudes  (28a)  as  (29a)  with  any  known 
solution  at  (12a)  for  e*(t),  and  q(()  =  Ce(rr  —  *)-  Analytical  expressions  for  c(c)  can  be  obtained  using  the 
local  analytical  solutions  (28b, c)  in  (28a)  or  (29a)  which  can  then  be  determined  by  the  method  of  stationary 


phase. 


V.3  Cross  Sections  from  Classical  Theory 

A  flsssicsl  theory  —in  the  sense  at  summation  over  probabilities  rather  than  of  amplitudes  so  that  inter¬ 
ference  effects  on  ignored—  is  readily  deduced  from  the  local  approximation  (26)  without  direct  evaluation 
of  c(«)  from  (28a).  The  theory  follows  from  (28)  with  conservation  of  probability  imposed  is.  from  the 
requirement  that, 

p4(t)  +  £  p{t)A  =  l  (30) 

where  P4  is  the  probability  je*|a  of  survival  on  V4  from  tx  —  0  to  t  and  where  P(()dt  is  the  probability  for 
autoionisation  within  the  time  interval  [t,t  +  A].  Thus, 


p(t)  =  =  ££!p„(t) 


(31) 


Within  this  time  interval  [t,  t+dt]  or  range  (ft,  ft+dft]  of  internudear  separations,  an  electron  is  emitted 
with  energy  in  the  range  [e,  t  +  de]  about  e(t)  =  W(t)  =  V1-K+  with  probability, 


P(t)  A  =  P(R)  dR  =  P(c)  dr  =  |c(«)|*  dr  (32) 

The  probability  densities  per  unit  interval  dR  or  dr  arc  then, 
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W-a-Sft-r(-ir5SH-S  «-> 

or, 

p(<) = woi*  *  -pw  ||p  s  p(x)  m*)rl  =  (ss*) 

respectively,  where  the  radial  apeed  *(R)  =  A,  W  =  dc/dt  and  W'(R)  =  dVF/diL  The  croaa  section  for 
disposal  of  all  possible  energies  <  in  A  -  B*  collisions  at  energy  E  is  then. 


<rj(E)  =  2r  J**  bdbj**  P(t)  dt 

=  2r  f*  b  A  i**  P(R)dR 
Jo  JMJt) 


when  the  mairiinum  impact  parameter  bx  which  just  penetratse  the  reaction  aone  ia  given  by, 


(34o) 

(346) 


£(*)  =  J$[l-2!iM]  (34c) 

This  croea  eaction  b  for  the  full  reaction  involving  all  <  and  all  ion  states,  ie.  for  both  Associative  and 
Penning  Ionisation  processes  Upon  Reintegration  (34b)  with  (33a)  reduces  to, 

which  is  the  standard  dassiral  result  for  any  absorptive  reaction.  For  small  autoionisation  widths  I*  which 
an  only  weakly  dependent  on  6  (or  L)  then, 


<m<) 

when  d*  is  an  element  of  the  trajectory  *  =  *{E,  6)  traced  in  time  dt  =  d*/u(R)  =  dR/v(R)  in  terms  of 
the  radial  and  local  speeds,  »(*)  =  A  and  u(R)  =  {2/U)l/*[E  -  V„(  A)]1/J.  respectively.  Since, 

-  2*  jf *  f*M  fW**  =  A*I*‘  **  l1  ~  m  dR  (34/) 

the  cross  section  (34e)  reduces  to, 

«<*>-?  (^rr^hspr*  «** 

in  agreement  with  the  results  of  Miller.  When  the  path  length  A,  —  hv(R)/T(R)  towards  autoionisation  is 
independent  of  R  then  (34d)  reduces  to, 
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'  'M-sr  *('-%)  {m) 

a  mult  reminiscent  of  Thomson's  tot  tarmolecular  ion-ion  recombination. 

Electron*  of  a  gives  energy  t  can  be  emitted  at  different  timee  t  j  and  tj  os  the  incoming  leg  and  again 
at  rx  —  i\  and  Tx  —  <*  on  the  outgoing  leg.  Thus,  W(t)  si  *  baa  four  roote  (aee  Figure  1)  for  i  in  tbe  range 
(0,  rjr).  The  total  probability  density  from  (33b)  tor  electron  emiaeion  is  then  the  aum  of  tbe  individual 
probabilitiea, 

pm = E  *«*»  *  t  S  r*  1  p*w  -  E  iwwr ‘  *(*)  (35«) 

where  the  total  probability  of  survival  on  the  incoming  and  outgoing  lap  at  R  is, 


i/'aaj 

*  /*.  *(*) 

Is 


(3tt) 

(35c) 

(35d) 


The  croas  section  tor  (A  -  B*)  collisions  for  disposal  (via  autoionisation  and  vibrational  excitation  of 
AB+)  of  energy  c  in  the  range  c ,  <  +  dc  is  <r(E,  c)  dt  where  the  differential  croas  section  is. 


~  *  **(*,«)  =  2*  jf**  *(«,*)»  *  (M) 

and  P(«,5)  is  given  by  (35)  which  implicitly  depends  on  the  impact  parameter  b  via  the  rlaaeical  Mbit 
R  —  R(E,  L)  where  the  relative  angular  momentum  l*  —  {2ME)&.  Not  only  is  t  the  vertical  energy 
separating  t^(JZ)  and  V+(N)  at  R,  but  it  is  also  (ct  |T)  the  vertical  separation  between  E  and  the  energy 
of  the  vibrational  level  *  of  AB*  (cf.  Figure  2).  The  probability  that  AB+  is  left  in  a  vibrational  level  v, 
assumed  to  form  a  quasi-continuum  is  therefore, 


J»(*.*).  *(«,») |£|  (37) 

Since  the  radial  action, 

J  =  fp+(R)dR  =  (*+0*  (38) 

is  quantised  for  bound  vibrational  motion  with  local  momenta  p+{R)  between  the  classical  turning  points, 


d*  _  dt_dJ_ 
dv  dJ  dv 


(39) 
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whan  v  ia  tin  frequency  for  vibrational  motion  in  level  *  within  the  potential  V+(R).  The  probability  for 
nutoioniention  with  AB+  left  in  level  •  ia  then, 

P(«,»)  *  4w*  g|VmW«))|*  A  (^y)  Pa(*)}  (40) 

which  ia  itimmainnleea  The  corresponding  croea  aection  for  eaaodntive  ioniantion  into  AB*(t)  ia  then, 

**(*.')  =  **/'*  PfaW*  (41) 

The  angular  momenta  L *  =  (2 ME)b?  —  (f + l/2)*h*  and  (J  +  l/2)*h*  of  relative  nuclear  motion  before 
and  after  autoioniaation  are  ao  large  in  compariaon  with  the  angular  momenta  of  the  ejected  electron  that 
t  —  /.  The  croea  aection  for  aaaociative  ioniantion  with  AB+  left  in  vibration al-rotatiocal  level  (v,  J)  ia 
•imply, 


=  rf-(2J  +  l)P(v,KJ))  (42) 

The  transition  (T)-matrix  element  ia  therefore  given  by  PJ^*(»,  b)  of  (40)  in  thia  rlaaairal  theory.  The  term  in 
bracce  in  (40)  is  a  -1-- 1  representation  (is.  excluding  phaaea)  of  the  Ftanck-Coodon  Factor  for  bound-free 
vibrational  transitions  (c£  Appendix). 

From  detailed  balance. 


«in*Va/(*;«.<0  =  (3w+)(W  +  l)kfoM(r,*,J)  ,43) 

where  and  w+  are  the  electronic  atatiatical  weights  of  AS*  and  AS *  and  2  is  the  apin-etatiatical  weight 
of  the  incident  electron,  the  croea  section  fore-  -  A2?+(v,  J)  dissociative  recombination  (ie.  for  the  reverse 
reaction  in  (13))  is  then, 


»«*(«.!.,  i)= ( jjf )  » { [iw(  *)f 1  (^)  p.(*)]}j 


(44o) 

(44k) 


where  the  total  survival  probability  P,{R)  on  the  incoming  and  outgoing  legs  at  R  is  given  by  (35b).  It 
will  be  shown  (cf.  |VI)  that  the  summation  in  (44b)  ia  a  representation  of  the  quanta!  matrix  element  for 
the  molecular  autoioniaation  frequency,  evaluated  by  the  method  of  stationary  phase  and  summed  over  the 
probabilities  (rather  than  the  amplitudes)  of  the  contributions  from  the  various  regions  of  stationary  phase 
at  R%.  The  term  in  braces  in  (44b)  is  the  contribution  at  to  the  Franck-Condon  overlap  between  the 
bound  and  continuum  vibrational  wavefunctions. 
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Moreover,  rnlmtioa  of  «,(<)  of  (29*)  by  tl be  method  of  stationary  phoao  yields  the  classical  path  T- 
matrix,  aop(«)  *  (h**)c,(t).  Neglect  of  interference  effects  between  tho  various  contributions  e*  to  tbs 
amplitude  e  from  each  region  •  of  station *ry  phase  yields  the  Hsssicsl  probability, 

P(e;h)  =  (**)5>|*.  (45) 

i 

given  by  (40).  Further  physical  insight  into  the  associative  ionisation/diseoristive  recombination  process 
(13)  and  simplified  expreseiona  for  the  respective  cross  sections  with  interference  effects  included  are  now 
obtained  by  evaluating  the  integral  (39a)  directly  by  the  stationary  phase  approximation. 

T.4  Stationary  Phase  Amplitudes 

The  main  contributions  to  the  integral  (39a)  arias  from  thoaa  time  intervals  surrounding  points  t,  of 
stationary  phase  given  by  y  =  dy(c;t)/dt  =  0.  The  phase  variation  outside  these  regions  of  stationary  phase 
is  then  sufficiently  rapid  to  justify  extension  of  both  integration  limits  in  (39a)  to  infinity.  For  one  region 
around fy  (say)  the  upper  and  lower  limits,  on  the  change  of  variable  toss  (t-fy)  becomes*  =  (t-ti)  — *  oo 
and  *(  =  — (fj  —  tx)  — *  —oo,  respectively.  The  integrals  in  (39a)  for  the  amplitudes  are  of  the  form, 

A*(s)s*  f(r,t)eap(±n(e;t))  dt  (46) 

For  cases  involving  at  most  two  stationary  points  tj  and  fy  which  correspond  to  phase  minima  (with 
y(r, fy)  >  0)  and  phase  maxima  (with  5(*ify)  <  0),  respectively,  insert  in  (46)  the  expansions, 

7(0  =  7(*i)  +  7 (t*)(t  -  fy)  +  )fr(W  -  <*)*  (47a) 

and, 


»(0  ■  »(fy)  +  *(fy)(<  -  fy)  (474) 

about  each  isolated  stationary  point  fy  given  by  y  s  dy(<;  t)/dt  =  0.  The  integral, 

£ **»  (*?w**)  d*=  [r]  '  ■*” (±*1)  m 

then  provides  the  Stationary  Phase  Evaluation  of  the  integrals  (46)  as  the  linear  expansions, 


A*(«)  =  Oi(«)  exp  [±<71  +  J)]  +  «»(«)  exp  [±<T*  -  J)]  (48o) 

or  equivalently, 

d*(«)  =  (oi(«)  T  «j(«)  «*p(±»7ai  )]  «*P  [±<71  +  J)[  (465) 
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of  the  individual  amplitude  with  magnitudes, 


•<(«)  *  [a»/  *(«).  »  =  1. 1 


and  phaae  differences. 


1*1  =  71-71  *?(*»)  “?(*»)  (48d) 

at  aach  phaae  point  4  treated  in  isolation. 

Note  that  the  relative  phaae  7*1  of  the  amplitude  O)  to  «i  in  the  term  A4  ie  reduced  by  w/2  to  give 
the  factor  of  — »  in  (41b).  The  phaae  change  from  1  -*  3  in  general,  ie  given  by  ^[aifivfi  -  sifwyi)*. 
Theee  expreeaiona  are  valid  for  well  aeparatad  regiona  of  Stationary  Phaae  ie.  for  771  >  1.  In  the  present 
application  t(« ,t)  ia  given  by  (2Tb)  tucb  that  ir(c,t)  s  (W(t)  -  «]/h  ia  the  only  derivative  of  7  that  depend* 
on  «.  Expressions  (48)  ate  therefore  appropriate  only  for  that  energy  which  satisfies  7  =  0.  When  =  W/k 
vaniahea,  ae  tut  an  extremum  value  e*  ia  W(t)  at  timet*  when  the  two  poiata  tu  of  atationary  phaae  coelaace 
at  t  (aa  for  a  cauatic  or  rainbow),  then  expression  (41c)  diverge*.  On  extending  the  eipanaimn  (47a)  for  7  to 
include  the  next  non- vanishing  term,  tf^t*)  <  0,  aince  J(t")  =  0,  and  on  recognising  that  t(c,  f)  =  (**-«)/* 
in  the  neighborhood  of  the  cauatic  at  r*  evaluation  of  (48)  than  yielda  the  finite  amplitude*. 


**(«)**  r~~r  ii<,r)Ai(rM)'mr±rr(rtr) 

*(*•) 


in  terms  of  the  Airy  function  At  of  argument, 


Expressions  (49a)  are  in  theory  valid  only  in  the  c-rangt  surrounding  the  caustic  at  c*  but,  in  practice,  only 
at  e*. 

The  analysis  above  is  formally  identical  to  the  well  established  analysis  of  daaeiral  rainbow  scattering; 
where  y(t),  W(t),t  and  c  above  are  analogous  to  the  phase  shift  7(f),  deflection  function  *(0  =  dt)/dl, 
angular  momentum  t  and  scattering  angle  #,  respectively,  in  elastic  scattering.  The  "transitional”  Airy 
approximation  (49)  does  not  uniformly  connect  with  the  "primitive"  result  (48).  By  mapping  the  phaae 
7(e;t)  onto  the  integrand  of  the  Airy  function,  a  Uniform  Airy  Approximation  which  uniformly  connects 
(48)  and  (49)  at  r*  ia  well  known  from  previous  work.  The  result  is  written  here  in  compact  form  as  the 
linear  combination, 


A+(e)  =  oi(<)ccp  [»(7i  +  J)]  F*(t«)  +  o»(*)e*p  [*(7*  “  J)]  F(t»i) 
■  («i(«)F*(7»)  -  •*»(«)  «*p(*7»i)F(tjx)J  **p  [*(7i  +  J)] 


and, 


A" («)  =  «i(<)  op  [-»(ti  +  J)]  r(ju)  +  •»(«)  sap  [-1(71  -  J)]  F*(tn)  (Wc) 

S  [•i(«)^(7>j)  +  *aj(<)  •»p(-»7ji)*‘*(7»i)]  **?  [~*(7»  +  j)]  (50d) 

when  the  complex  function  F  in  defined  in  terma  of  the  Airy  function  Ai(i)  and  itn  r -derivative  Ai'(s)  fay, 

FIt*i(«)]*  [*i/V'4Ai(-n)  +  «r*/>*-^Ai'(-n)]eV-»(2S.- J)  ;  j |x|*/l  =  7.1  >  0  (51) 

for  Tai  =  7*  —  7i  <  0.  Since  711  in  the  area  encloeed  by  the  W(t)  carve  and  the  straight  line  W  a  «  (cf. 
Figure  1)  it  in  always  positive,  except  when  it  in  sera  at  c  =  c*  =  It  in  shown  below  (|T)  that  the 

divergence  in  the  constructive  interference  tern  (•>  +  *1)  at  the  caustic  a  0)  is  exactly  balanced  by  the 
vanishing  of  the  coefficient  s1^4  of  At ;  also  the  divergence  in  the  coefficient  m~1**  of  Ai*(-s)  at  the  caustic 
is  offset  by  the  destructive  interference  term  (at  —  aj)  which  vanishes  mote  rapidly.  In  the  limit  of  high 
s  >  1,  or  for  well  separated  regions  yjt  >  1,  F’(tji)  — »  1,  with  unit  amplitude  and  aero  phase  such  that  (60) 
tends  to  the  primitive  form  (41).  The  Uniform  Airy  result  (SO)  is  general  in  that  it  continuously  connects 
the  Caustic  (tu  *  0,-fj  =  0)  result  (40)  at  <*  with  the  result  (40)  for  writ  separated  regions. 

Autoion  i motion  Amplitude  and  Interference  Patterns 

Application  of  the  above  Stationary  Phase  Prescription  (50)  to  the  integrals  (30a)  for  the  autoionisation 
amplitnds  c(<)  yields  after  some  raarraagment, 


te(t)  a  ['Pj/,(€)PJ1  +  *Pj/,(r)#y1  sap(-«7n)  +  Vi/a(t)Fit  enp(-m»)  +  «*rf-n«)]  *•?  [-*(71  +  J)] 

(52°) 

which  form  illustrate!  how  the  prngrsasive  contributions  are  added  as  time  cvoules.  An  equivalent  form  is, 

*K«)  =  [^i/J(«)  {«t  *n  +  «*?(-»74i)}  +  vPa1/J(f)  «ep(-»7ai)  {cj/Vj  -  tsj/ji  «*p(-ma)}]  «*P  [-*(71  +  j)] 

(525) 

which  illustrates  the  addition  of  the  inward-outward  contributions  at  each  JL  The  probability  densities  in 
(52b)  for  autoionisation  at  the  isolated  points  Ri  are, 

*(0=x|v*W))|3*|l^|  1  =  ^^|wr'(^)l"1-  -<«**.*>  (**) 
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In  (52*),  th«  probabilities, 


Vi  *  V>  =  P*<4(X*) 


=  a);  ^4  =  Pi«2(*i) 

already  incorporates  the  amplitudes, 


(52d) 


and, 

for  survival  from  Jly  to  it  on  the  incoming  leg  or  from  Rx  —  R<  —  Roa  the  outgoing  leg,  respectively.  The 
indices  »  =  1  —  AmVi  are  associated  with  quantities  calculated  at  times  ti,ta,ts  =  ry-(j  and  t|  *  rx  - 1| 
which  correspond  to  Si  ^  lij  "*  H)  Hi  on  the  incoming  and  outgoing  kgs  of  the  rleeeical  path  as  t 
progresses.  The  phase  differences  in  (52a)  for  times  on  the  same  leg  are. 


7><  =  jfjTVW-*1  *  (52») 

and  for  times  tj  >  U  on  different  legs  are, 

w-  jjT*  *+zj TVw-«]  *  (*»> 

where  t,  is  the  time  of  closest  approach.  Thus  (52)  is  the  classics!  remit  (35)  fully  generalised  by  this 
classical  path  theory  to  incorporate  relative  phases  The  autoionisation  probability  density  P(<)  =  |c,(t)|J 
will  therefore  exhibit  various  interference  effects  arising  from  the  phase  differences  between  the  amplitudes 
contributing  at  the  same  Rt  on  the  incoming-outgoing  legs  and  at  the  different  R\  and  X]  on  either  the 
incoming  and/or  outgoing  leg*  For  well  separated  regions  7x3  >  1  and  (52a)  tends  to, 


c(«)  =  esp(— *731)  +  Pj/J(c)  s*p(-n»x)  +  »'Pi/,(r)  e*p(-i74i)  (53o) 

where  the  (redundant)  multiplicative  phase  factor  +  J)]  of  (52a)  has  been  neglected  on  the 

RHS.  The  full  interference  pattern  exhibited  by  the  probability  density  is  then, 

P{t)  =  |e(«)|*  =  +  V9  +  V»  +  VA  +  2(P1?j),/*  «» TSi  -  2(*>,P,)1/»  sin  7*3 

+2(PiVi)lf*  tin 7«  +  2(P»7>i)l,i  liny*  +  2(ViPt)l,icosjn  -f  2{P%P^)x,i  eo$y^  (54) 

=  |  J*x /8(«)  {«!  +  **i  **p(— *74i)}  +  *P*1,a(«)  {e*  «ap(-*T*i)  -  IS*  «»p(-mi)}( 


This  illustrate*  nrioui  iatdmaei  tfbeU  arising  (a)  from  diitwat  Si  mad  J!i  oo  the  bum  lag, 

Plt  =  Vi  +  V,  +  2(PiP,)l/*  *11.7*1  (55a) 

P*  =  Vt  +  V 4  +  2(W)l/*  (55k) 

(b)  from  the  in-out  contribution*  at  a  given  S», 

Pi*  =  V1  +  Vi  +  2(7>,7>4)l/*  tinni  (Wo) 

ft*  S  P,  +  Vi  -  2(7>i  V»)l,*nnyu  (55k) 

and  (c)  from  different  J?j  and  JSa  on  different  lege, 

P«  »  +  v*  +  2(PiP*)1/*  00*7*1  (57o) 

Pu  =  Vt  +  V t  +  2 (PaP*)l/*  oo*7 4a  (57k) 


All  of  the  above  interference  pattern*  Py(c)  oscillate  with  c  about  the  classical  mean  fP<(e)  +  P>(«)) 
between  the  envelopee  (f>11/*(«)  ±  P,l/,(f)]*  with  frequency  2whji(t).  Sine*  7<i(«I  k)  »  *l*o  a  function  of  the 
impact  k  or  angular  »"«■"—*«"«  £  S  /,  the  interference  pattern*  will  be  effectively  waahad  out 

in  the  k-integration  for  cross  sections  each  a*  Penning  loniiation  in  (24).  For  croa*  aectiona  aa  associative 
ionisation  or  dissociative  recombination  involving  AB+  with  specified  rotational  state  /  ( ie.  at  a  selected  im¬ 
pact  parameter  k)  the  above  interference  patterns  will  persist.  The  Uniform  expansion  for  the  autoionisation 
probability  density,  provided  hr  the  present  theory  (52*)  is, 

P(e)  =  J  [P11/,(«)P«  +  <Pt*(*)Pii  «*p(-*7Ji)  +  i  **p(-*m)  +  »Pi/,(«)FJi  esp(-i74i)]  |(5«<s) 

=  |  [p,l,*(e)  {eiFai  +  i*i  JJi  ««p(-*74i)}  +  *P}%)  { ca^ai  **p(-*ysi)  -  **»*at  «*p(-mi)}]  j  (55k) 

Here  in  addition  to  the  (rapid)  -oscillations  which  vary  with  e  and  k,  broader  oscillations  due  to 
the  Airy  functions  within  Pji  will  appear  in  the  vicinity  of  the  caustic  at  e  =  «*.  These  oscillations  are 
independent  of  k. 


▼.6  Classical  Path  Cross  Sections  for  Ionisation/Rocombination 
la  tbs  reaction, 


A*  +B++AB*  +•-  (59) 

tbs  energy  «  =  V4  —  V+  of  disposal  is  distributed  between  tbs  energy  t,  of  tbs  ejected  electron  sad  tbs 
energy  **t  abeorbed  by  tbs  molecular  km  AB*  in  vibrational-rotational  level  (**/}.  The  square  of  tbs 
angular  momentum  of  tbe  A*  —  B  orbital  motion,  with  impact  parameter  5  is  l%  as  (Slfxa  =  /(/+1)A*, 
and  is  conserved  throughout.  Tbe  integral  cross  section  for  associative  ionisation  foe  disposal  of  all  «  is 
therefore. 


*,(*)  =  2*  jf  *  *  *(•!*,»)*  («0) 

where  JP  =  |o(*)|J  is  given  by  (SB),  and  where  ij,  tbs  maximum  impact  parameter  which  results  in  tbe 
distance  Rx  of  orbital  closest  approach  is  given  by, 

(«) 

Tbe  (-integration  is  over  all  (singb  valued)  energies  <  between  0  and  tbs  maximum  possible  u*  Bar  a  specified 
separation  Rm  =  R(* m)  (cf.  fig.  1).  The  cross  section  (60)  for  population  of  all  accessible  (n,f)  states  of 
AB*  including  the  vibrational  continuum  is.  for  both  Penning  and  Associative  Ionisation  is  then. 


"W-jiahj i.  U'L 

(«2a) 

s-a- r~ (2i  +  l)di  r~  P(^;n,f)dn 

"AM  Jtm 0  JO 

(625) 

s5?;S>*+1)£  TOM) 

(62c) 

» 


Tbe  probability  for  ejection  of  an  electron  with  energy  c,  =  E  ~  tnt  and  with  tbe  ion  AB *  left  in 
vibrational-rotational  level  (n,f)  is, 


P(J;n,f)  =  P((;*rI)i 
where  tbe  vibrational  energy  spacing  (39)  is, 


(«3) 


d*  _  d«  dJ 
dn~  dJ  in 


(M) 
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in  terms  of  the  quantised  actios  J  *  fpitdR  =  (n  +  l/2)h.  Hence, 

P(f;n,f)  =  |e(c«,«){a  (*»w)  =  W«t  £«•)!*  («) 

with  —  E  —  t •*.  The  amplitude  or  the  classical  path  transition  matrix  (Tc^|*  obtained  from  (U)  ie 
therefore, 

*>(»)  =  IW 

=  K»)|*  =  +  <P,1/a(»)*Ji  «*fK-7ii) + 


(M-) 

=  |P*/a(»)  [ciFu  +  leiJPJj  esp(— net)]  +  »PJ1/*(»)  «*p(-*7ii)  [«»**»  -  o«p(“nn)]|*  («*) 

The  probabilities  and  Pi  are  related  by  (52d)  and, 

«<")  =  (*««(<.)  =  {ma)|-  (^y)}-  (W) 

ie  dimensionless.  The  croae  section  for  autoionisation  and  population  of  AB+  is  level  (a,/)  is  therefore, 

*ju(X;n,l)  =  jr-(2f+  1)  W»»,«».)l*  *  (M) 

where  a(n,f;c«)  defines  the  classical  path  T-matrix  Tap.  From  the  detailed  balance  relation  (43)  the  cross 
section  for  dissociative  recombination,  the  inverse  of  (89)  is, 

*»*(«•)  =  (;j£+)  jr  W*i  *  *«)l*  (•») 

Since, 

^  |  V.  (*(.))!•  (7.) 

then  on  defining  the  Hx— “dependent  partial  cross  section, 

•«-bIs(20  FPHwfo)}  <”•> 

for  capture  at  the  stationary-phase  point  R,  the  cross  section  can  now  be  expressed  as, 

»nx(«5 »)  *  J»J/a  lejfti  +  ««iFJi  esp(-vy«i)]  +  «*p(-*Tai)  («**>!  -  **riat  «*p(-mj)]|  (Tib) 
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when  s  v(Ri),  *  =~ i e(Hi)  and  as  =  ^(Ri)  el  (52e,f).  This  is  the  besic  erpwssion  for  the  cross  section 
for  dissociative  recombination  in  the  present  dsssicsl  peth  theory.  For  one  point  of  stationary  phnse  then 
Fji  =  1  end  9%  as  0  in  (71b)  which  redaces  to, 


***(*;  n)  =  a{R)  (cj(*)  +  wJ(JI)  -  2c*(*)a,(B)  rin y«i]  (71c) 

end  which  exhibits  the  expected  in-oat  interference  st  it  In  (71b)  end  in  the  key  results  (S2b),  (58b)  end 
(68b)  above  the  phase-differences  given  formally  by, 

7 a  =  l£ {W  -  [v+(t)  +  e)}  *  (71  d) 

when  Ri  =  R(ti,U)  end  Rt  =  £(t],ts)  an  determined  in  practice  as  follows.  Since, 

iu4(R)  =  =  S-  (r+(H) + «]  -  ~  (72) 

and 


i  V.J(  *)  =  ^jlSi  .  *  -  K,(«)  -  fg-  (7J) 

where  the  nuclear  angular  momentum  J  ~  L  remains  conserved,  then, 

W(X) -«  =  V4(R) -  |V+(S)  +  <] 

=  £[*♦(*) -*2(*)I 

■MM*)  -*<(*)]» 

when  an  averaged  (common  trajectory)  radial  speed  is, 

V*<  +  *+)  dR 

2  if  “  dt 

The  stationary  phase  condition,  «  =  W(Jt),  therefon  demands  the  conservation  of  kinetic  energy  of 
the  nuclei  (ie.  k4(Ri)  —  k+(Ri)).  Both  the  angular  momentum  (L,J)  and  kinetic  energy  of  nuclear  relative 
motion  an  therefon  conserved  in  the  CP-theory.  In  terms  of  the  phase, 

A(*)=  f*[k+(R)-kt(R))  dR  (75) 

Jk. 

measured  relative  to  that  at  R*,  the  distance  of  closest  approach,  the  phase  j, 

y(t)=i£m*)-']  *  (76) 


(74o) 
•  (745) 
(74c) 

(74 d) 


is  expressed  as, 


fA(Jlx)- 

\a(*X)  + 


-  &{R)  t  <  t. 


The  phass  differences  is  ell  the  previous 


A(H)  t  >  t, 
n a  (62b),  (S«b)  and  (66b)  are 


calculated  in 


Tji  =  Ai  -  A*  s  Am 
Tm  *  Ai  +  A} 

7M  =  2Aj 
741  =  2Ai 


when  A<  =  A(JZ<).  The  claaaical  path  amplitude, 


oox(»)  =  <k')1/,c(f) 


when  c(c)  ia  given  fay  (52a)  ia  therefore, 


«CJ»(n)  =  {^i/a(»)  («ifti  +  «i*ii  e«p(-2»Ai)J  +  «P,l/,(n)eap(-iAu) 

-  *^fti eap(-2»A,)]} cap  [i(Aj  -  J)]  esp(-iA (Hx))  (766) 

when  J\(«)  an  given  by  (67).  Thia  expreanon  (76b)  ia  alao  confirmed  by  a  stationary-phase  evaluation  of 
(29a),  aince  the  pfaaaaa  (29b)  and  (75)  an  identical. 

For  one  region  of  ataiionary  phaae  (78b)  ndncea  to, 

«cx(»)  *  2 wVZ(X)  {|1T(JI)|  ( jg-)  }1/l  MR)  +  MR)  e*p(-2iA(*))]  (79) 

which  ia  in  agreement  with  the  aemiclaaaical  reault  of  Miller. 

The  prcaent  claaaical  path  theory  (71b)  and  (78b)  therefore  furnish ea  quite  naturally  the  uniform  gen¬ 
eralisation  to  two  regions  of  stationary  phase  including  the  caustic  region. 

Special  Cate t:  Turning  Point  and  Caustic 

Two  cases  of  special  interest  arise.  One  is  when  Ri  =  X(<i)  and  Rt  =  R(U)  coalesce  at  the  turning 
point  Rt  where  the  radial  speed  v(R)  vanishes  and  the  cross  section  (71a)  diverges.  The  other  case  is  when 
Ri  and  Rj  coalesce  at  the  caustic  wherein  c*  =  W(R)  it  «n»Timnm  (cf.  Figure  1)  so  that  vanishes 

and  (71a)  again  diverges. 


Turning  Feint  Diverge nee:  The  divergence  it  J 2«  cu  he  tnUid  by  ideptiaf  la  (71a)  the  ipwd, 


in  the  presence  of  autoioaisation,  rather  than  -*  0  u  with  Ifilkr  Q.  A  peasant  proposal  for  dealing  with 
tba  divergence  at  tba  turaiag  point  k  baaad  on  tba  recognition  that  an  average  over  tba  quantal  distribution, 

K(*)f  ** *  =  M  -  .  (SI) 

can  be  replaced  by  the  corresponding  rlsssiral  average  over  tba  period  T*t  for  vibrational  motion  in  level 
n  =  v  so  that  2*w/*j»  in  (71a,c)  is  simply  f^+(Jl)|*.  Then  (71c)  with  cufJZ*)  =  s*(A«)  and  7*1  =  0,  yields 
tbs  finite  result, 

"*('•1 - [isb^rKSf)]  {t'Wwwfarj-ijf  ™«*])  m 

»«r.(*.)ft(*..*x)  (828) 

This  naturally  decomposes  into  a  croaa  section  r,  for  capture  at  JZ«  and  tbs  probability  P*  toe  survival 
to  tbe  point  Rx  of  stabilisation.  Ibis  simple  result  then  represents  a  first  improvement  over  that  of  Bardaley 
(cf.  equation  (4))  in  that  it  includes  W'(R)  =  [V*  -  V+]  ratber  that  VJ  alone.  The  original  result  (4)  is 

therefore  valid  (a)  when  one  region  of  stationary  phase  at  J2t  is  assumed,  (b)  when  V*  is  so  shallow  that 
&+(£«)  =  0  and  (c)  when  V4  is  so  steep  that  tba  Winans-Stueckelberg  wavefunction  |VJ(.R<)|-1^*  f(R  -  R,) 
can  be  used  (see  Appendix)  for  the  continuum  vibrational  state. 

The  divergence  at  JZ«  in  general  originates  from  tbe  normalisation  of  tbe  JWKB  wavefanctions  implicit 
in  the  vibrational  overlap  S  in  the  absence  of  autoionisation.  The  term  in  braces  in  (71a)  is  a  JWKB  ap¬ 
proximation  to  this  overlap.  This  divergence  is  eliminated  when  Airy  functions  rather  than  JWKB  functions 
are  used  —  even  when  the  effect  of  autoionisation  on  the  normalisation  is  ignored.  See  the  Appendix  for 
further  discussion. 


Caustic;  The  divergence  in  (71a)  or  in  (67a)  due  to  the  sero  in  W'(R)  at  the  caustic  JZ(ti)  =  Rite)  is  exactly 
balanced  by  the  behavior  of  tbe  function  P(tji)  in  (66a)  as  tjj  -*  0.  This  becomes  apparent  by  rewriting 
the  contribution, 


«(«)  =  J^Wnm)  +  U»,1/a(n)F*(7»»)e.p(-mi)  (83) 

to  the  amplitude  (66a)  from  the  incoming  leg  in  the  alternate  form, 

•<«)  *  [p*/a(»)  +  **'*(»)]  A(7,i)  +  .  [P,l/*(n)  -  Pj/a(«)j  A'(tm)  (64) 


32 


where. 


A(t*i)  *  w*/*e1/4A*(-e);  *  =  [W4|*/a  (85) 

A'(t*i)  =  »1/*e-l/4Ai'(-e)  (M) 

The  equivalent  forme  (S3)  and  (84)  for  a  ere  ueeful  for  probing  the  separate  limit*  Tn  >  1  end  Jti  -*0 
,  respectively.  In  the  neighborhood  of  the  caustic  at  i*,  where  <*  =  W(t*)  is  maximuTn  and  W(?)  =  0  then, 


so  that  c  s  W(t)  at  the  two  times, 


ta,i«l*±  [2(«*-  «)/|^(f)|]1 
The  derivative  W(t)  therefore  tends  to  sero  at  <i,j  ae, 

[*(«*-‘)  |^(**)|]Xi 

The  phase  is  expanded  consistent  with  (87)  as, 


t « = Jin + amt  -  n + \  (^r)t  (<  -  n*  (*>®) 

*  t(<-) + £(«*  -  «)(<  -  r)  -  gf  |*(f)|  (t  -  r )•  '  (*») 

since  y(t)  =  (W(t)  —  e]/A  and  since  $(i*)  vanishes.  The  phase  difference  at  lj,i  is  therefore, 

2»/t  fe*  -  *\*/» 

1».M  «t(«.)-7(i.)-  7iT|»(,.)| 

The  argument  of  the  Airy  functions  in  (84)  is, 

..  |  2ft  \l,$  («*-«) 

The  (c*  —  e)1/4— dependence  of  m1^4  in  (85)  is  exactly  balanced  in  (84)  by  the  (<*  —  c)~ll4— dependence  of 
1  in  (67)  of  Pi1/*(s).  Also  £p*^*(*)  -  P,1/J(c)]  tends  to  sero  foster  with  t  than  The  classical 

path  cross  section  (71b)  for  dissociative  recombination  then  reduces  to, 


*D*(e)  =  <r(ff)  [<£(*)  +  «3(P)  -  2cs(ff)*s(JI)  sw2A(ff)] 


where, 


33 


(•3*) 


■  •£:(*!*)  l1?1]  {t  la^jl 

which  is  finite  at  the  caustic.  The  teem  in  braces  represents  the  dominant  contribution  from  the  caustic  to 
the  Franck-Condon  Factor  |5t*  -in  effect  replacing  the  corresponding  (divergent)  term  in  (71a). 


T.8  Correspondence  between  Classical  Path,  Semiclaseical  and  Quantal  Thaoriae  of 
Dissociative  Recombination 


In  order  to  probe  with  the  quantal  result,  correspondence  the  symmetrical  term  (29a)  of  the  banc 
classical  path  result  is  useful.  On  making  the  substitution  (74),  the  rlsneics  1  path  amplitude  a<7p(n)  = 
where  c  is  given  by  (29a)  is  therefore, 


ecj»(n)  =  2» cap  [— tA(Rx)]  J  V£(.R)  |(v/h)l/*  9~l  (e*(.R)  cap  +«A(R)  +  «(*)  cap  -iA(R)]}  dft  (94) 

The  term  in  braces  is  (94)  is  ssesntislly  the  product  of  the  bound  and  continuum  vibrational  wavufunctions 
(eee  below).  Stationary-Phase  Evaluation  of  (94)  yields  (78b),  as  it  should. 

Quanta!;  The  quantal  exprsaeion  for  the  autoionisatioo  frequency,  with  electron  energy  s  in  the  range  c,  c+dc 
and  with  the  ion  left  in  state  (w,  J)  is, 

I  n*,{r'R)  *  '«*>)!«)*  (*> 

where  the  system  wavefunctioa  for  A  -  B*  collisions  at  energy  E  is, 

Si  =  R)Y, *(&)+&,  R)  (98) 

the  product  of  the  Boro- Oppcahsi mar  electronic  wavefunctioa  ds  the  actual  continuum  (radial)  vibrational 
wavefunctioa  de  in  the  pteaence  of  autoionisation  and  the  rotational  wavsfunction  Y>*r  (A).  The  rovibrational 
wavefunctioa  for  AB+(w,  J)  is. 


*+{&)  =  ~4+(R)YJlt(A)  (97) 

Both  fi(E,  R)  and  the  continuum  electronic  function  ++  for  the  (c~  -  AB+)  system  are  energy  nor¬ 
malised  with  unit  densities  p(E),  p(t)  of  states,  respectively.  The  incident  current  ( 4j/dE)dE  integrated  for 
all  directions  of  £  is  therefore  (8 wME/h*)dE  =  (h^f/t^hidE.  The  associative  ionisation  cross  section 
{dt/Jdj)  is  then, 


°ai{E)  =  rf-(2/+  l)Pq(E, «;/,.) 


(88) 
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whan  the  quatal  autoionisation  amplitude  or  transition  matrix  element  ia, 

•«(*)  *  2*  jT  »£(*)  [tfW-Wj  ** 

By  detailed  balance,  the  dissociative  recombination  czoee  section  ia, 


(M) 


(100) 


which  is  dimensionless. 


Rtdnciio*  of  Quontol  to  Stmieiotneol:  The  product  of  the  eemidassical  JWKB  vibrational  wavefunctiona, 


for  the  diacrete  levels  with  inner  (left-hand)  tuning  point  Bo  end, 


for  the  continuum  lsvds  subject  to  autoionisaiion  is. 


(101) 


w*’  *  j^JSjp  (£*'"+ 1) dm) 


It* (*hM*)  *  (*W*)1/,(*+»4)*l/*  M&)  e«p+tA(Il)  +  *o(B)c*p-iA(B)] 

where  the  phase  dHEweare  is. 


(103) 


A(«)=  /  k+(R)dR-  fMk4[R)iR  (104) 

•  r,  «/*. 

The  highly  oeriUatocy  exponentials  of  the  phase  sums  have  been  neglected  since  they  provide  little  relative 
contribution  to  (99).  The  quanta!  amplitude  (99)  with  the  semiclaesical  (103)  is  then, 

eq(e)s2w(^/h)1/>jf* (»+»r)-1/,V2,(B)  [ct[R)  oxp  +wh(*)  +  s*(*)  top  -»A(K)]  dR  (105) 

The  stationary  phase  condition  yields  h+(Jt)  =  hs(Jt)  so  that  the  arithmetic  mean  f(iZ)  and  geometric 
mean  [»+('R)»r(f*)]t^s  are  identical  at  the  points  of  stationary  phase  and  (105)  is  identical  to  the 
path  amplitude.  Stationary  phase  evaluation  of  but  (105)  and  (94)  confirms  that, 


oq(»)  =  ocp(n)  o*p+tA(Rx)  (106) 

where  «cp  is  given  by  (71b).  The  rlassical  path  amplitude  (94)  is  therefore  the  quanta!  amplitude  (99)  with 
the  JWKB  vibrational  product  (103),  thereby  establishing  the  equivalence  between  the  matrix  element  (105) 
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ud  it*  Fourier  Componsat  RtpnmUtiot  (Nt)  or  (29a).  Tbii  equivalence  for  vibrational  transition*  is 
•omswhai  okia  to  tbs  Heisenberg  CorreeponeUxxc*  Principle  for  electronic  transitions 

Wboa  tbs  bound  vibrational  JWKB  warsfunction  is  tskw  with  rsspsct  to  its  right-hand  turning 
point  Ho  thon  tbs  analysis  is  os  abort  but  with  A  replaced  by, 

A(H)  *  £  *+(*) «  -  jT  *<(*)  **  (107) 

which  rsducas  with  tbs  aid  of, 

k+(R)dR  =  (•+ 
to 


Tbs  smaller  of  A  and  A  in  practice  would  bo  adopted  in  tbs  general  result  (79b). 

7.7  Rats  of  Direct  Dissociative  Recombination 

For  a  Maxwellian  distribution  of  electron  energies  s  =  <*(97)  at  temperature  T,  tbs  DR-rate  is, 

o(T)  =  ¥  r  oDk(*Y  «*PH'  &  3  (vdkCT))  (1W) 

where  ¥  is  the  mean  electron  speed  (SAr/vifan)1^*  and  odr  is  a  mean  crons  section  at  temperature  T. 
An  energy  threshold  «o  s  I^(JZq)  —  V+(R ,)  >  0  pertains  for  the  case  when  the  energy  V^Ay)  at  the 
crossing  exceeds  the  original  vibrational  energy  of  AB+(v).  In  terms  of  the  probability  (M)  for  dissociative 
recombination, 


A(R)  =  (v  +  i)»-  J^*+(R)dH  +  jT*MJl)« 


(106) 


«(T)  =  (;■£*)  jT  J^njO^-r/hTJdr  (110) 

is  the  bade  expression  for  the  rate  which  includes  all  the  phase  information  in  (66).  In  order  to  obtain 
a  simplified  analytical  rate,  assume  that  there  is  only  one  region  of  stationary  phase  at  Jt,  so  that  c  = 
V, i(R')  -  V+(Rt)  and  that  the  in-out  interference  effect  can  be  ignored.  The  probability  is  then, 


w«««)P = W(«.) +w 


(ilia) 


which  with  (35)  is, 


w** = **•*{«■*  [-±£*  fgf «] — [j  £  }  (>>■*) 
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when  tha  Urb  in  bnM  is  the  »vm(i  probability  for  mninl  from  S«  21  Rj  tad  li,  ^  fi,  Ij  on 
V4(R).  Adoption  of  (111)  ia  (110)  therefore  demands  the  location  R*  which  dapaada  ia  turn  on  tha  energy, 

*  =  +  (ill) 

of  relative  motion  uadar  V4.  From  a  =  V r4(R*)  -  V+{R«)  and  tha  vibrational-  rotational  land  (e,  L)  at  AB+, 
k+(Rf)  —  k4(R*)  aad  F  can  than  ba  datarmiaad  to  provide  £,ho  fuactioo  at  *.  Furthar  reduction  at  (110) 
with  (111b)  ia  therefore  not  poaaibla  without  additional  assumptions. 

(a)  la  tha  low  aaargy  limit «  -*  0,  (£,,  R,)  -*  Rx  and, 

Wnja  -  0)|*  «  2*r(Rz)  |^(Hx)f‘  {  [^y]  ■  fc+(*x)|*}  (US) 

auch  that  (110)  raducaa  at  low  T  to, 

*(r>  “  (aJa-W*  {l*(«x)|,*r|»’(jfe)|-}  <»«> 

Tha  tarm  in  bracaa  ia  an  effective  Franch-Condoo  factor  for  bound-froa  transitions.  Thia  analytical  raault 
ia  a  generalised  version  of  tha  original  raault  of  Bataa  ia  that  it  include  W(R)  =  —  V+)  rather  than 

V't  alone.  It  therefore  allows  for  a  distinction  to  ba  made  between  craaajnga  on  either  aide  of  tha  potential 
minimum.  For  croaringa  at  the  potential  minimum  tha  results  are  identical. 

(b)  Assume  either  that  V*  ia  so  shallow  that  1 +(£«)  ~  0  or  that  Vs  ia  so  steep.  Than  R,  =  R,  which 
is  given  universally  by  a  as  V4(R,)  —  l r+(Rg).  Tha  recombination  at  electrons  at  aaargy  a  origins  tee  at  the 
distance  of  closest  approach  ao  that  tha  cross  section  (82a)  can  be  used  directly  ia  (100)  to  give, 

where  P(c)  ia  the  probability  density  (33b)  given  by, 

where  c  =  W(R).  fat  constant  drainage  dP4/dt  is  constant  R,  and  Rx  and  (114)  is  recovered. 
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App«ndbu  Franck- Condon  Overlap  for  Bound- IVoe  Vibrational  IVimMmm  with  and  without 
Autoioniaation 

1.  JWKB  Franck- Condon  Overlap  without  Autoionisetion 

The  JWKB  normalised  aamidaaaical  w* refunction*  for  the  bound  vibrational  level  (n,  f)  of  AB+  with 
vibrational  frequency  v  ia, 

=  (Ai) 

where  h*>  =  ia  the  level  spacing,  and  Bo  ia  the  rlearirel  turning  point  given  by  the  ianermoat  aero  of 

3*<(*)  =  =  *-  (*+(*)  +  «)-  ( A2 ) 

the  radial  apeed  v+(B)  of  relative  motion  of  energy  (f-i)n  potential  V"f(B).  The  JWKB  wavefunctioa 
energy  normaliaed  to  S(B  —  F)  to r  the  vibrational  continuum  of  AB*  without  autoioniaation  ia, 

*<(*)  *  Pm^t^***  [jC  *<(*)*Jl+ ?]  •  *»*«  (") 

where  fi,  ia  determined  by  the  innermoat  aero  of, 

^  »  *  -  W)  -  (jm) 

for  the  radial  apeed  «e(B)  of  relative  motion  ia  the  dissociative  potential  V*(B).  Angular' momentum  of 
relative  nuclear  motion  ia  conaerved  (J  =  X).  The  Franck-Condon  amplitude. 


SwjTd;(  B)da(B)dB  (A5) 

ia  then  written  with  (Al)  and  (AS)  aa, 

5  =  (x) 1,3  £  [«*p+»A(B)  +  enp -*A(B)]  dB  (A8) 

where  the  phase. 


A(B)  =  /\(B)dB-  [*  kt(R)dR 

Jr*  Jr. 


haa  a  stationary  point  where  A'(B)  =  dA/dB  =  0  ie.  where  k+(R*)  —  4*(B«).  The  previous  condition 
(equation  (38)), 


Va(B,)wV+(B.)-M  (A8) 
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for  4  vertical  transition  at  root  R  m  R,  m  then  recovered.  The  (±)  terms  **p(±iA)  is  (At)  provide  the 
contributions  to  S  from  tba  (— )  and  wt|oia|  (+)  componnia  of  d 4 .  On  expanding, 

A(*)  *  A(B.)  +  A'(*«)  (*-*«)  +  jA"(Jl.)  (*  -  R,)*  (At) 

and  on  changing  the  integration  variable  to  a  =  R  -  R,  with  limits  (±00)  then  (At)  can  be  evaluated  with 
the  rid  of, 

j~~  **P  (±*  M  •*)  do  =  [j^J  esp  (it J)  .  ( A10) 

The  Stationary-Phase  Vibrational  Overlap  is  thsreforc, 

-  w  l/'^)[isw]  *  a  ♦  -  -  h*>  *  a }  «“> 

where  the  constant  phases  (±*/4)  pertain  to  positive  or  negative  values  of  A "(H.)  =  k\(R,)  -  h'd{R,), 
respectively  ie.  to  either  minima  or  mamma  in  A  at  JS«.  Item  (A2)  and  (A4), 

a"<*>  *  s^bs  f*  <AU> 

where  W  =  V*(.R)  —  V+(Jl).  The  Frenck-Condoo  Factor  is  then, 

W  *  [^y]  mR,)fl  mn *  [AW  =F  J]  (AW) 

which  oscillates  (rapidly  for  A  large)  about  its  average  value, 

^**yb],iFwri  (xw) 

For  one  root  12,  of  (At)  whan  W  s  dW/AR  <  0,  the  (+)  sign  in  (AW)  is  appropriate.  Several  variations 
of  (AW)  can  be  constructed. 

(a)  In  classical  mechanics  tbs  quantal  probability, 

I^Wl’dJlrc^avdt  (A15) 

is  replaced  by  the  corresponding  classical  average  over  the  period  T  for  vibrational  motion,  the  factor  of  2 
arising  from  inward  and  outward  radial  motion  ie.  |^*(J2)|*  =  2v/o(R).  This  also  follows  from  the  JWKB 
function  (Al).  Use  of  this  correspondence  in  (Alt)  therefore  yields  the  Fraack-Condon  Factor, 

\S\\  *  2 1<(£,)|*  |^(Ve  -  V +)  1  rin>  [a(*.)  =F  J]  (Alt) 

which  is  advantageous  in  that  it  dreumvents  the  divergence  in  the  overlap  (A13)  at  the  classical  turning 
point  R,  common  to  all  JWKB- based  approximations.  Then 
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|5l!  =  K(*.)|*|^(V4-v+)|^ 

which  k  finite. 

(b)  The  simplest  energy  normalised  continuum  wavefunction  it  the  Winans-Stueckelberg  function, 

\dVd\~U* 

MR)  =  1^1  t(X-X*)  (A17) 

when  fi,  ia  the  classics!  turning  point.  The  Franck-Condon  overtop  if  then, 

W*«WW|*|2r|  1  (A18) 

I  lit* 

to  be  compered  with  the  more  accurote  ixpr— ion  (A16)  or  (A14).  It  k  therefore  valid  when,  \%t\  ^  j^gf1 I 
ie.  the  potential  V4  k  eo  steep  and  strongly  repulsive  relative  to  V+(R).  Tbs  Frenck-Condon  overlap  (A1S) 
k  that  used  fay  Bardaley  to  provide  the  cross  section  (4)  for  dissociative  recombination. 

(c)  Airy  Function  Remedy 

In  order  to  remedy  the  well-known  breakdown  of  the  JWKB  functions  (Al)  and  (AS)  close  to  the 
classical  turning  points,  the  JWKB  functions  can  be  replaced  by  their  Airy  function  counterparts, 

em^j£i<«+^  |s»'*«jf^idR  (A19) 

in  (Al)  and  (AS).  The  stationary  phase  result  (A1S)  k  than  replaced  by, 

\s\*  =  y^y]  I  r(jyr‘  [»*' V'*Ai(-n)]  *  (A20) 

for  the  resulting  overlap,  where  the  argument  of  the  Airy  function  A*  in  terms  of  the  phase  difference  (A7) 
k, 


^)s:[|A(R.)]a/*. 

For  large  arguments  t),  ie.  for  R,  well  removed  from  the  classical  turning  points  JZ<  and  JZq,  then, 


(A31) 


*lfW/4AH-Tj)'-^  m  (a  +  J)  R<>Rc  (A22) 

so  that  (A1S)  k  recovered  (for  the  case  W'(R«)  <  0).  The  overlap  (A20)  uniformly  connects  the  classics! 
■■■saaible  and  inaccessible  regions  and  does  not  diverge  when  R,  k  located  at  the  classical  turning  point  JZ«. 

2.  Franck- Condon  Overlap  with  Autoloniaation 

As  noted  in  f  XXC,  the  effect  of  autoionisation  on  the  A  —  B *  system  k  realised  by  sitifning 

the  complex  potential  energy  B, 
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(AM) 


*<(*)  =  v„(*)-jr(ji) 

to  the  decaying  dkoditin  state.  The  JWKB  solution  of, 

=  [*- Vd(R)  +  Jr(Ji)]  s  SfcllhMJl)  (A24) 

which  satisfies  the  appropriate  conditions  for  the  incoming  end  outgoing  wane  at  the  boundary  Rj  is, 


W(^)  =  re  /  D\ll 

[ku(R)Y 


("s) 


where, 


ij(K)  =  *J(R)  +  »^r(JZ)  (A2So) 

W  <  *  -  Vs(*)  (X2M) 

The  autoionisation  is  in  effect  within  the  reaction  soae  between  the  crossing  point  Rx,  where  Vd  =  V*, 
and  the  distance  R.  of  dosast  approach  at  energy  B.  The  physics  within  da  abort  becomes  apparent  by 
rewriting  (A25)  as, 

I**"  ("ijT  $>  ■"')*•’’  [-  (£““+ f)] 

-^[-s(r+£)™"h*K*'"+j)i  ,"7*) 

*  (ESP  5)  -  *">■*•  (£*•"+ 5)]  ("7‘> 

where, 

is  the  survival  amplitude  from  JZx  to  /Z  on  the  incoming  leg  and, 

<x2,) 

is  the  amplitude  for  in-out  survival  for  the  sequence  JZx  -♦  R*  JZi-  The  stationary-phase  determination 

of, 
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5  =  y  <•(*)*-(*)**  (A30) 

with  two  poiata  of  wtll  separated  lUtionvjr  phase  when  h+( A,)  =  h^A*)  at  Ai,j  yields, 

S  =  (Si  {ci  +  tat  e*p(-2iAj)}  +  tSt  esp(-tAu)  {e*  -  *•*  eap  (-2iA»)}]eapt(Ai  -  (AS la) 

when, 

^  *  [k5I)  |wr*(Ji‘)rl]l/*  {A5lh) 

when  Ci  =  <*(&),  *  =  *e(Hi)  end, 

S  A(J2i)  —  A(Jfy ) s  A,- Ay  (ASlc) 

in  terms  of  A  given  by  (A7).  The  notation  h  each  that  Jlj  >  J?j  when  W'(Ri)  <  0  and  W'(.Rj)  >  0, 
respectively.  When  autoionisation  is  neglected,  c*  —  1  =  ae,  then  (A31)  reduces  to  (A13). 

For  two  general  regions  of  stationary  phase  then  the  Stationary-Phase  Uniform  Prescription  (SO)  yields, 

S~  [Sl{cIFi1  +  ieiFJ1e*p(-2*Ai)}  +  i5ieap(-»Ala){eJFJ1-ie^,aie*p(-2«AJ)}]enpi(Ai-j)  (AS2) 
for  the  Franck-Condon  overlap,  and  where, 

/■[Ai^  =  [ir1/,sl/4Ai(-s)  +  *sJ/,s“1/4A»f(-*)]  np-t  -  j)  ;  j  |e|*/#  *  Au  >  0  (ASS) 

is  the  function  introduced  previously  in  (SI).  The  above  results  (A31)  and  (A32)  for  the  vibrational  overlap, 
with  and  without  (c*  ss  1  ss  as)  autoionisation  for  two  general  regions  of  stationary  phase  appear  new  to  the 
literature,  and  help  establish  the  correspondence  between  the  present  clawical  path  theory  and  the  quanta! 
result. 
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n.  THE  MODIFIED  ADIABATIC  INVARIANCE  METHOD 
FOR  THERMAL  ION-DIPOLE  MOLECULE  REACTIONS 

If.  R.  Flannery  aad  X.  Qi 
School  of  Physics 
Georgia  Institute  of  Technology 
Atlanta,  Georgia  30332-0430 

Abstract 

The  somirlasskal  adiabatic  invariance  method  (AIM)  is  used  to  compute  thermal  energy  rate  coefficients 
for  capture  of  various  atomic  aad  molecular  ions  by  diatomic  and  triatomic  molecules  with  permanent  dipole 
moments.  Particular  ion-molecule  systems  studied  in  this  paper  include  Ee+,  C+,  aad  E%+  ions  reacting 
with  polar  SCI  molecules;  and  H*,  C+,  0+,  ECO*,  and  Be*  ions  with  HCN  molecules.  In  addition, 
a  new  modification  of  the  adiabatic  invariance  method  is  presented  in  this  paper,  which  accounts  for  the 
coupling  of  the  internal  rotational  angular  momentum  (j)  of  the  target  molecule  to  the  orbital  angular 
momentum  (f)  of  the  projectile  ion  about  the  center  of  mam  of  the  target  molecule.  Comparison  of  the  AIM 
results  with  the  present  modifications  to  the  adiabatic  invariance  method,  aad  with  available  experimental 
data  aad  extent  theory,  for  the  above  ion-molecule  reactions,  indicate  that  the  inclusion  of  j  —  t  coupling 
within  the  AIM  is  moat  important  in  the  limit  of  light  ion  mam  or  weak  ion-molecule  interactions. 


PACS  no.:  34.50.Lf,  S2.30.Fi,  82.32.Nr 
Running  Title:  Thermal  Ion-molecule  Reactions 

Preferred  Section  in  Joumek  Molecular  Interactions  and  Reactions,  Scattering,  Photochemistry. 
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VI.  1  Introduction 


The  ser nklasrical  adiabatic  invariance  method  (ADI)  vu  ftnt  proposed^1)  for  ion- molecule  collisions 
when  the  Interaction  potential  energy  ie  an  adiabatic  invariant  about  the  rotational  angle  of  the  target 
molecule,  and  haa  been  applied  to  ion-molecule  colliiione  with  linear^*)  and  eymmetric  topW  moleculee  with 
pennenant  dipole  momenta  and  ion-quadrupole  molecule^4!  colliiione.  However,  at  prevent  no  account  haa 
been  made  in  the  above  AIM  calculation!  for  the  coupling  between  the  internal  orbital  angular  momentum, 
j  of  the  target  molecule,  and  the  relative  angular  momentum  f  of  the  projectile  about  the  center-of-mass  of 
the  projectile  ion-target  molecule  system.  In  thia  paper  a  modification  of  the  invariance  method 

ii  propoeed  which  taken  into  account  the  coupling  between  j  and  t  in  ion-molecule  ceBSaom 

An  outline  of  the  remainder  of  the  paper  ie  aa  follows.  In  |ZX.A.l  the  modifications  to  the  ADI  are 
presented  (hereafter  denoted  AIM(j,i)),  while  in  XX. A. 2  an  overview  of  the  original  AIM  ie  given  in  the 
notation  of  thia  paper.  The  detaila  of  the  computation  a t  the  capture  cruae  sections  and  rate  coefficient!  into 
specific  and  quant*!  states  of  the  target  molecule  are  given  in  |XX.S.  A  discussion  of  the  present  results 
with  available  theory  and  experimental  data  is  provided  in  1XXX,  while  a  summary  of  tbs  paper  ie  given  in 
{XT.  Unless  otherwise  noted,  all  quantities  are  in  atomic  units. 

VI. 2  Theory 

2.A.1  Modified  Adiabatic  Invariance  Method 
The  Hamiltonian, H,  describing  the  reactions 

« 

X*  +  AB  —  {XAB+Y  -  products  (la) 

X*  +  ABC  -*  (XABC*  )*  -  produets  (lk) 

of  an  ion  X*  with  diatomic  (ALB)  or  triatomic  ( ABC)  target  molecules  with  permanent  dipole  moment  D, 
in  the  body-fixed  CM  frame,  ie 


h-b?  t  S±  |  Dco*9 

;  2pr>  2r4  + 


(2) 


where  B  and  04  are  the  rotational  constant  and  dipole  polarisabilhy  of  the  molecule  respectively,  and  with 
r,  9  representing  the  intermolecular  distance  and  angle  between  the  dipole  moment  D  and  the  intermolecular 
axis  of  AB,  respectively.  The  interaction  potential  for  the  ion-molecule  reactions  (1)  is, 


*M)— fl  +  Sfj r* 

Since  /  is  a  conserved  quantity  in  the  collision  (1),  we  can  write  (f-  j )*  in 
(C)  of  /  onto  the  intermolecular  axis, 


(3) 

of  the  projection^ 
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(4) 


{f-  J)*  *  (J  +  j)*  +  J*  -  *0* 

The  Hamiltonian  H  can  then  b«  written, 

,<  =  j’J,  +  5^<-'+S>,+?-s0’]-|s+£5:1-  <6> 

The  adiabatic  potential  is  defined  for  fixed  r  by 

*W  *  (B  +  2^a)j*  +  («) 

where  psp*  +  pv*/#m*d.  Since  p  is  cyclic,  p,  ss  m,  the  integration  of  p  over  a  complete  cycle  (range) 
of  angle  8  is  an  adiabatic  invariant, 

+  (7) 

To  determine  ti«  adiabatic  potential  c(r)  from  the  quantisation  condition  (7)  on  pe,  an  mn— inn 
relating  p»  to  e  is  needed.  Such  an  exp r— ion  is  obtained  from  the  definition  of  j*  (given  above)  in  ter—  of 
Pi  and  py  and  making  nee  of  the  fact  that  p  is  a  cyclic  variable  and  results  in. 


"•  =  (b+5^s)  *[*> 5^)^»]  '  (,) 

Substituting  (8)  into  (7),  and  defining  the  dimensionless  variables,  p  =  cos#  and  s  =  ry/WfD,  and  with 
« (*)  =  t/B,  equation  (7)  can  be  rewritten  as, 


w  '  '  r  (2pD**  +  i)t  J,  '  i — p* 

where  we  have  replaced  n  on  the  left-hand  side  of  (7)  with  /—  |  m  |  due  to  the  invariance  of  the  integral  of 
Pi-  Note  that  when  the  limit  r  — *  oo  is  taken,  the  integral  in  (9)  can  be  done  analytically  with  the  result, 
n  +  1/2  =  v/’H0O)~  |  m  |=  (j  +  1/2)-  |  m  |.  The  function  /  in  the  integrand  of  equation  (9)  is  given  by, 


/(*)  -*V*  -p  + 


((*- »v-^) 


=  (p  -  o)(p  -  h){p  -  c) 

where  the  roots  {a,  k,e}  of  the  cubic  polynomial  are  ordered  such  that  a  >  b  >  e.  After  so—  algebra,  the 
integral  equation  (9)  can  be  rewritten  as  a  transcendental  equation  for  the  dimension!—  adiabatic  potential 
«/"*(*). 
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U-  W  +  j)» 


.*+ 


1 

2mX> 


(«“*)*  =  -  «)X(t)  +  2(a  -  e)JF(f) 


(11) 


when  JT,  X,  and  II  represent  the  complete  elliptic  integral  of  ftnt,  second,  and  third  kind^,  respectively, 
with  moduli  defined  as,  4  =  (5  -  c)/(a  -  c),  pi  a  (h  -  c)/(l  -  c)  and  p»  *  (e  -  h)/(l  +  c). 

The  adiabatic  potential  ujm(s)  ie  then  obtained  from  (11)  once  the  roota  of  the  cubic  polynomial  (10) 
are  determined.  In  addition,  to  determine  the  capture  croee  section  in  fIZ  .1,  an  expeesrioa  ie  required  tot 
the  first  derivative  of  u,m(s)  with  respect  to  s.  By  applying  Leibnits's  rule  to  (9),  the  required  expeeaeion 


for  the  firet  derivative  ie  obtained, 


da  ~ 


2a"* 


(• 


(— (*+,> 


(12) 


2.1.2  Adiabatic  Invariance  Method 


In  this  section  we  provide  a  short  overview  of  the  original  adiabatic  invariance  method!1-*!  using  the 
present  notation  for  completeness  In  collisions  of  ions  with  molecules  with  permanent  dipole  momenta,  the 
adiabatic  potential  is  defined  fin  find  r  by. 


<r)  =  Bj1  + 


Deoti 

“P— 


(IS) 


The  problem  is  again  to  determine  the  adiabatic  potential  using  the  quantisation  condition  (7).  Ne¬ 
glecting  the  coupling  between  j  and  t,  the  analog  of  (9)  is, 


«0-W+s)  = v7(d  TT7  <“) 

with  the  function  /  in  the  integrand  of  (14)  given  by, 

/(*)  =  A*  -  a*!!^  -  P+  («-  m*)»* 

(IS) 

=  (p  -  o)(p  -  *)(p  -  c) 

where  the  roota  are  ordered  a  >  b  >  e,  and  the  remaining  terms  in  (14)  and  (15)  are  the  same  as  in  (9)  and 
(10).  Equation  (14)  can  be  converted  from  an  integral  equation  to  a  transcendental  equation  in  w  with  the 
result, 


U  -  H  +  “  «)*  *  2(«#m(a)a*  -  a )K(q)  +  2(«  -  e)S(q)  - 


m*n* 


n(pi,t)  n(pi,y) 

1-c  1  +  c 


(19) 
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where,  in  the  cm  vhaa  K  <,B  the  captor*  Croat  atetioa  oim  k  replaced  with  the  second-order  Stark  croaa 
section 

= -  (f )  <») 

when  the  term*  Aj m  an  equal  to  a  -  Cjm  if  a  >  Cym  and  aero  otherwise.  The  conataata  C>m  above  an 
equal  to  -1/S  if  j  =  0  and  /  if  j  >  0,  while  a  =  atB/D2. 

The  rate  coefficient  for  capton  into  the  atate  of  the  molecule  k  then  given  by,- 

Upon  averaging  over  an  aaramed  MaxweU-Boltsmann  thermal  distribution,  an  averaced  rate  coefficient 
for  capture  at  tamperatun  T  k, 


£f  (3j  +  l)h*(T)eap-*U+i/t>V» 

(MT)>  =  i»  -  (M) 

f(Ji  +  l)aqp“*0+»/*)*/w 
f-o 

when  the  sum*  over  j  range  from  taro  to  jmsm  which  k  the  maximum  orbital  angular  momentum  of  the 
target  molecule  which  atiil  aupporta  an  angular  momentum  barrier  in  the  affective  potential. 

VI. 3  Results  and  Discussion 

The  molecular  parametsn  of  HCl  far  as,  D  and  B  used  in  the  praaent  calculations  .wen  3.93  (A*), 
1.09  (Delya),and  10.99  (cm-1),  respectively,  while  far  ECN  the  corresponding  values  used  wan  3.59  (A*), 
3.99  (Helps),  and  1.49  (cm-1).  The  adiabatic  potential  energy  for  the  present  modified  AIM  results  wen 
determined  from  aquation  (11),  while  the  original  AIM  results  shown  wen  obtained  from  aquation  (19). 
The  complete  elliptic  integrak  of  the  first  and  second  kind  were  computed  from  standard  codas,  while  the 
complete  elliptic  integral  of  the  third  kind  was  computed  by  the  method  given  by  Carlson  W.  The  capture 
cross  sections  ey(JJ)  were  calculated  from  (33)  over  an  energy  range  B  from  0  to  3000JB.  The  averaged  rate 
coefficients  far  capture  wen  calculated  from  equation  (35)  with  jm—  equal  to  30  for  ECl  and  40  for  ECN. 
The  weighted  error  in  the  rate  coefficient  for  HCl  for  jmm  =  30  at  1000  *K  k  less  than  0.3  percent,  while 
the  error  for  ECN  far  j — r  =  40  at  1000  *K  k  leas  than  3%. 

Table  1  gfves  the  comparison  of  the  rate  coefficients  of  the  AIM(j,f)  and  the  AIM  to  experimental 
data  and  other  theories  for  capture  for  H*,  C*,  O*,  and  ECO *  ions  reacting  with  ECN  molecule  at 
temperature  305, 300, 440,  and  540  *JT.  Table  3  and  3  present  the  rate  coefficients  for  Bf  ion  reacting  with 
ECN  and  ECl  molecule  at  temperature  of  305,  300,  440,  and  540  *JT.  The  experimental  data  shown  in 
tables  1-3  were  taken  from  reference  [9].  Table  4  shows  the  rate  ^efficients  of  Be*  and  C*  ions  reacting 
with  ECl  molecules  at  37,  99,  and  300  *K  with  the  experimental  data  taken  from  reference  [9]. 
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For  the  Hi ction  E* +HCN,  the  rate  coefficients  computed  from  the  AHi(j,/)  and  AIM  ut  significantly 
different  (aw  figure  1).  It  indicates  that  the  rotational  sad  orbital  angular  momentnm  coupling  ( j—t  coupling 
)  Ass  effect  during  the  recombi nation  process.  His  AJM(j/)  results  an  close  to  the  results  of 

the  sdisbstk  capture  and  centrifugal  sudden  approximation  (ACCSA)J10*  The  results  of  both  ADi(j/)  and 
the  ACCSA  agree  well  with  experimental  data  except  at  205  *JT.  The  cause  of  the  overestima tattoo  by  the 
theories  is  still  not  dear. 

The  rate  coefficients  of  reaction  of  B$  ion  with  ECN  and  SCI  molecules  are  well  predicted  by  the 
AIM(j/)  and  AIM  as  shown  in  figures  2  and  2,  respectively.  The  small  difference  between  the  results  of  the 
AIM(j|f)  and  AIM  indicates  that  the  j  —  t  coupling  does  not  play  an  important  role  during  the  collision.  In 
'Cables  S  and  4  excellent  agreement  of  the  AIM(j,f)  and  AIM  results  can  be  seen  with  the  clesssical  trajectory 
method  (CT)  statistical  method  •  probability  weighted  effective  potential  method  (PWEPM),  flux 

weighted  effective  potential  method  (PWXPM)  l11-14-1*!,  quantum  mechanical  method  •  statistical  adiabatic 
channel  model  (SACM)  and  ACCSA  at  temperatures  above  20  *1. 

For  the  reaction  ffe+  4-  MCI,  table  4  and  figure  4  show  that  there  is  virtually  no  difference  between 
the  rate  coefficients  computed  by  the  and  the  ADI.  The  j  —  l  coupling  in  this  case  is  entirely 

negligible  .  thereby  confirming  that  the  j  —  t  coupling  is  important  only  for  vary  fight  ion  or  weak  ion* 
molecule  interactions.  The  results  of  both  the  AJM(j/)  and  the  AIM  show  excellent  agreement  with  the 
ACCSA  and  the  SACM  results.  All  theories  agres  with  the  measured  data  to  within  the  experimental  error. 

For  the  reaction  of  C*  with  ECN  and  MCI  molecules,  ths  AIM  results,  along  with  other  theories,  over 
estimate  the  rate  coefficients  when  compared  with  experimental  data.  Studies  by  Clary  at  a L,  indicate 

that  the  overestimation  is  due  to  the  open  shell  nature  of  the  carbon  ion  C*  in  the  collision.  With  at- nitric 
calculations,  they  suggested  that  for  C*  ths  final  results  should  be  multiplied  fay  a  factor  of  2/3.  Hence,  in 
this  paper,  the  present  modified  AIM  nsults  for  the  case  of  C*  have  been  multiplied  by  the  factor  of  2/3. 

A  shwiliar  situation  exists  in  the  case  of  the  reaction  0+  4-  ECN.  The  overestimation  of  extant  theories 
in  the  case  of  reactions  of  oxygen  ions  0+  may  again  be  due  to  the  open  shell  nature  of  0+  as  was  indicated 
fay  at- initio  calculations  for  the  case  of  C+  ions.  At  present,  no  comparable  ek-***tie  calculations  have  been 
done  for  the  case  of  oxygen  ions  reacting  with  ECN. 

The  rate  coefficients  for  the  reaction  ECO +  4-  ECN  computed  by  ths  present  AIM  calculations  agree 
with  the  experimental  data  and  the  ACCSA.  At  205  and  300  *ff  the  praaent  AIM  results  are  about  25 
percent  above  the  experimental  data,  which  is  within  the  experimental  error  ban  of  the  SIFT  technique. 


VX.4  Summary 

A  new  modification  to  the  adiabatic  invariance  method  (AIM(j,f))  has  been  presented  and  shown  to  be 
important  for  thermal  ion- molecule  reactions.  The  principal  feature  of  the  modification  to  the  AIM  theory 
is  ths  inclusion  of  ths  coupling  between  the  internal  rotation  angular  momentum  of  the  target  molecule  with 
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the  orbital  angular  momentum  of  the  projectile  ion  about  the  CM  of  the  target  molecule 
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Thbb  I.  Comparison  of  rata  coefficients  for  uptun  for  tlu  mctloa  X*  +ECN  (X  =  H,CtO,MCO) 


Ibapcratnre  (*K) 

205 

300 

440 

540 

reaction 

method* 

k  (xlO***^#-1) 

B*  +  HCN 

AIM  04) 

(present  work) 

IS. 

IS. 

10. 

M 

AIM 

IT. 

15. 

12. 

11. 

(present  work) 
ACCSA 

16. 

13. 

11. 

11. 

experiment 

11. 

11. 

10. 

6.7 

C*  +  HCN* 

AIM 

S.9 

3J 

2.6 

2.6 

(present  work) 
ACCSA 

S.T 

3.1 

2.7 

2i 

“*"ri“-‘ 

S.4 

3.1 

3.0 

2.9 

0+  -f  ECN 

ADI 

5.4 

4.6 

M 

Si 

(present  work) 
ACCSA 

5.0 

4.2 

3.6 

3.4 

experiment 

3.7 

3i 

3.3 

3.1 

ECO*  +  ECN 

AIM 

4.6 

3.6 

3.2 

3.0 

(present  work) 
ACCSA 

4.2 

3.5 

3.1 

2.9 

3.7 

3.1 

2.9 

2.6 

(a)  for  definitions  at  tbs  symbols  ns  table  2.  Tbs  experiment  and  ACCSA  data  were  from  [•]. 
(t)  the  AIM  and  ACCSA  remits  bare  all  beta  multiplied  by  2/S.  [IS] 
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Table  So  CwnpiAoa  of  rate  coefficient!  for  uptun  for  tko  B»+  +  HCN  reaction 

Temperature  (*JT) 


method 

15 

50 

100 

150 

205 

220 

280  280  300 

k  (10-,em**- 

320 

*) 

380 

440  500 

540 

600 

AIM(i4)a 

35.5 

10.7 

14.0 

11.5 

9.9 

9.8 

8.8 

8.5 

8.0 

7.4 

8.9  8.5 

8.2 

6.9 

AIM* 

38.8 

20.2 

14.4 

11.9 

10.3 

10.0 

9.3 

9.0 

8.4 

74 

74  74 

8.7 

6.5 

CT*  (A) 

31.8 

20.5 

8.7 

(B) 

22.0 

13.0 

74 

PTRAT1 

32.8 

18.8 

13.6 

11.3 

9.9 

9.8 

8.9 

8.7 

8.2 

7.8 

7.2 

8.6 

CVTST* 

53.5 

29.4 

20.9 

174 

14.8 

14.3 

13.2 

12.8 

12.4 

12.0 

11.1 

10.4 

9.4 

ADO* 

8.5 

5.4 

PWEPM* 

28.1 

14.4 

10.2 

8.5 

74 

74 

8.8 

8.8 

8.2 

5.8 

5.5 

5.1 

PWEPM* 

31.0 

17.0 

11.9 

9.9 

8.8 

8.4 

7J 

7.8 

74 

8.7 

8.3 

6.8 

ACCSA* 

64.0 

24.0 

12.0 

9.8 

We4r 

7.3 

6.9 

SACM* 

18.0 

11.0 

1 

7.2 

64 

opt.* 

9i 

8.1 

7.1 

8.7 

Prutnt  Work- 

(*)  AIM(j,4):  Modified  Adiabatic  laniiuet  Method 

(b)  AIM:  Adiabatic  Invariance  Method 
Clwieal  Tktorits: 

(c)  CT:  Claancal  T^jeetocy  (A)  [11]  and  (B)  [13] 

(d)  PTRAT:  Parametrised  Trajectory  [12],  [13] 

(a)  CVTST:  Canonical  Variational  Transition  State  Theory  [11] 

SUtUtieal  Tkatiu: 

(0  ADO:  Average  Dipole  Orientation  [•],  [18],  [17] 

(t)  PWEPM:  Probability  Weighted  Effective  Potential  Method  [11],  [14],  [15] 

(h)  PWEPM:  Flux  Weighted  Effective  Potential  Method  [11],  [14],  [15] 
Quantum  Tkuoriu: 

(i)  ACCS  A:  Adiabatic  Capture  and  Centrifugal  Sudd  vs  Approximation  [8] 

(j)  SACM:  Statistical  Adiabatic  Channel  Model  [18],  [17] 

(k)  ezpt:  Experiment  [8] 


Tfcblt  3.  CompariiM  ofrat*  cotffcfanti  for  captvra  far  tha  f»+  +  SCI  raaction 

Tkmpcraiun  (*JT) 


method^ 

15 

SO 

100 

150 

305 

330 

260  360  300 

k  (10-,cm*s- 

330 

1 ) 

380 

440  500 

540 

600 

AXMU4) 

13.3 

7.4 

5.6 

4.9 

4.3 

44 

4.0 

3.9 

3.6 

3.7 

34 

3.4  3.3 

3.2 

3.1 

ADI 

ISA 

7.6 

5.7 

4.9 

4.4 

4.3 

4.0 

4.0 

3.9 

34 

3.6 

34 . 34 

3.3 

34 

CT* 

13.8 

6.6 

5.4 

4.4 

3.9 

4.0 

3.9 

4.0 

4.0 

34 

3.4 

34 

34 

PTRAT 

13.6 

7.6 

5.6 

5.0 

4.4 

4.3 

4.1 

4.0 

3.9 

34 

3.6 

34 

34 

CVTST 

19.4 

10.9 

8.0 

6.6 

6.0 

5.9 

54 

5.4 

U 

54 

4.9 

4.7 

4.4 

ADO 

3.4 

PWEPlf 

9.9 

5.9 

4.5 

4.0 

3.6 

3.5 

3.4 

3.3 

3.3 

34 

3.1 

3.0 

2.9 

FWEPlf 

11.7 

6.6 

5.1 

4.4 

4.0 

3.9 

3.7 

3.7 

3.6 

3.5 

3.4 

3.3 

3.1 

ACCSA 

16.0 

7.6 

4.5 

4.1 

3.6 

34 

34 

SACK 

13.9 

7.6 

5.7 

expt. 

4.2 

3.9 

34 

3.4 

(t)  for  dsBiritions  and  w&imw  of  the  symbols  sss  tabls  3. 
(t)  CT:  Classical  Trajectory  [11] 
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TabU  4.  Caapttboa  of  rat*  cooficUnts  for  captur*  tor  th*  rwetioa  X*  +  ECl  (X  =  Ht,C) 


Ikmperaiure  (*JT) 

37  M  300 


reaction 

method* 

k  (10_,em*«_1) 

H*+  +  ECl 

AIM  04) 

(preeeat  work) 

8.7 

5.7 

3.4 

AIM 

(preeeat  work) 

8.0 

5J 

3.4 

ACCSA 

8.4 

5.5 

33 

SACM 

8.8 

5.8 

3.0 

11.  ±  3.0 

4.8  ±  1.4 

3.3  ±  0.5 

c+  +  jrci* 

AIM 

(preeeat  work) 

3.7 

3.5 

1.5 

ACCSA 

3.5 

3.3 

1.8 

SACM 

3.8 

3.3 

1.3 

experiment 

3.8  ±  1.1 

1.8  ±  0.8 

1.0  ±  OX 

(a)  for  definition*  of  th*  ajnnbob  aaa  table  3.  Th*  experiment,  ACCSA,  and  SACK  data  van  from  [0]. 
(f)  th*  AIM,  ACCSA,  and  SACK  reeulta  have  all  b*ea  multiplied  by  3/3.  [IS] 
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